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ABSTRACT 


The singular value Bode plot of return difference and loop 
gain matrices have emerged as useful indicators of 
multivariable robustness. The H, and H, control theories 
provide a systematic procedure for shaping the singular value 
loop gains of a multivariable feedback control system. 

It is shown that H, control theory, using specified 
performance objectives and stability constraints, is effective 
in synthesizing a stabilizing controller for the statically 
unstable longitudinal dynamics of the X-29. H, Conia 
synthesis also demonstrates a good ability to cope with a true 
multivariable design problem such as the multiple, 
independently controlled surfaces of a super-maneuverable 
aircraft. However, it is also shown that the control surface 
deflections and control rates necessary to effect the 
specified performance levels exceed the performance 
capabilities of the X-29's actuators. 

A work-around to the limited actuator performance is 
provided by penalizing the control input vector more heavily 
during the problem formulation. This approach, while reducing 
the actuator performance requirements, results in a limited- 


performance X-29. 


bts 


TABLE OF CONTENTS 


ieee TRODUCTTON Sec ce ct cc twee ec ee we ec ee et eet ee eects 5 << 1 
il. FEEDBACK PROPERTIES OF MULTIVARIABLE SYSTEMS ....... 4 
Ae LURN Sn PP ERENCE MAUR TX ieee o 6 oe TS we te eee wee + 
B. MATRIX NORMS AND SINGULAR VALUES ..........22.. 6 

C. MULTIVARIABLE ROBUSTNESS AND SINGULAR VALUE 
MO Cori A VUNG Beiitais 66.6 0 6 6 Sad eicles = 0 6 8 ce ee mc ee ce ewe 8 
Meme, CONTROL DESIGN 2.2... ccc were cece sc ccc cece crc ce ecnn 20 
Pee OD er Me PORMULATUON ee ew cw wt ee eee we ew eee ON 
Bow SELECTION OF WEIGHTS ........ Dueis ociciicie «os efeusmenele « < 28 
See ONrROL SYNTHESIS .......2. UMMM Solel ss. ee oMetene ls % 0 0 


D. COMPARISON OF LQG, H,, AND H, CONTROL 


PUN MEM NcIs sce 6c 6 cet tc hw eh te ee ee ewe ee 8 a2 

Pe eeormCcONTROGEER SYNTHESIS .....c.cccc ccc cece n es cee 40 
ee erm OUR MM MoMmoeREPTTON 2... eee ce wt ee ce eee es 41 
Fe tte) ee ©) EC eres < cee. - « ceMeetencie cs cls cet tec ee ee 43 

eee Gime Se NTEE STS (5. ss oo ic wc cc te cc tee wees 49 

Pee eoreGiernnoUbho OF THE H, CONTROLLER ....scce--- 62 
eth TED-PERFORMANCE H. CONTROLLER DESIGN ...... 84 

Pe niCrUSotONS AND RECOMMENDATIONS ........cccccccccs IOS 
Sees )|6 6A 29 MATLAB SCRIPT FILES ......-ccccccccces 108 
Seen ExX Bs X-29 MODEL STATE SPACE REALIZATIONS ....... 134 
MES Te iE > au amernct scaler. s «oc ccs 060s 0 ss ee eee eee eee 730 
hy BESO S10) SUM CIC). a ie SY eee 2 


LIST OF TABLES 
ORDERED LIST OF THE UNCOMPENSATED X-29 MODEL 
STATES «2. ce cccvcscceces erererrcemnrrrrmrrnr so 
UNCOMPENSATED X-29 OPEN LOOP POLES .........5.5.506.66 600 
H, COMPENSATED X-29 CLOSED LOOP POLES. 2. nee 


CLOSED LOOP POLES OF THE LIMITED=PERPOCRMANGES =. 


V1 


LIST OF FIGURES 


Pamela tense meomiellUuStraClOMes sess. ssa ve ec ee eee 5 
BEeCabDaCk, COMmerGimeoy Stem 2°. Seger s e e  e e wee. LO 


Properties of the Relative Error Bounding 
FP TWUAGIES OWES 6 5 6 6 6 I 14 


MIMO Singular Value Loop Shaping 


Requirements ........... Tee oc Ce 16 
Relationship Between o,,,(I + G(s)F(s)] and 

IMEC NCNM Ses eet ee ee ee 19 
Sanda rdeH. Sia lmiemGal mmr TOD Lem ammis vis6 cs sess esc ee 20 
Pee Chine 1) CCl Wiis Memes crs chara. SMEG es oe ee ee RIS. 24 
Singular Value Bode Plot of S, T, W,', 

and W; | IPN, GUNUsis less SG so + + liable so eel 
Meee atinie (PROCESS (oi... . os cw ewe se eee wee eae e aes Se 
PO Cmmee GMa YS COM ss << <6 600 0 © eilbeie © 0 SleliSiic Sicielte ss see oe B35 
ey nee Te cee CE a) meet. als 2 os ss ee ee 38 
Uncompensated X-29 Open Loop Configuration ......... 42 
Singular Value Plot of the X-29 Plant G(s) ......... 45 
Open Loop Bode Plot of X-29 83rd Order 

ieee eT MMIC s ec cc se sees cc es occ eet tt et te ee 8 47 
eee ee DeSIGiy SPECILECACIONS 6. fale ss cc 5s 3 5 5 We we ee 48 
Pere DCS tm opeemstacatlOnsS 2. . 5. . siiiie . - + silieieles . ss 49 
Pepeost Function ||amiy for yl ......62...2.05..- 51 
Eeecost Function |/T,,l|, for y=6.7 ........----.-.- 51 
Bepeost Function |/T. lla for y=12.5 ...........45.. 52 
Senotcivity Hunetion S(s) for H, Solution, y=1...... 53 


Sensitivity Function S(s) for H, Solution, y=6.7 ... 53 


Vil 


aalea 


ele 


wks 


~14 


eS 


wo 


al Fe 


~LS 


lo 


aU 


oz 


ee 


25 


~24 


Ges, 


26 


ey 


Zo 


eo 


Sensitivity Function S(s) for H, Solution, 


Complementary Sensitivity Function T(s) 


for H, Solution, y=1 . 


Complementary Sensitivity 
for H, Solution, 7=6.7 


Complementary Sensitivity 
for EO Solution, y=l2eo 


Comparison of 8th Order and 10th Order 
X-29 Controllers 


8th Order Controller Additive Error E, 


Comparison of 10th Order and 15th Order 
X-29 Controllers 


10th Order Controller Additive Error E, .. 


Feedback Configuration for H, Compensated 
X= 29 


Singular Value Plot I+G(s), Uncompensated 
Meo 


Singular Value Plot I+G(s)F(s), 
H, Compensated X-29 


Singular Value Plot I+(G(s)F(s))7', 
H, Compensated X-29 


Universal Gain and Phase Margin Curve..... 


Singular Value Plot I+F(s)G(s), 
H, Compensated X-29 


Singular Value Plot I+(F(s)G(s))', 
H, Compensated X-29 


Advanced Control Modes 


Compensated X-29 a and gq Responses to 
TPs 


Compensated X-29 a and q Responses to 
inet 


Compensated X-29 Control Deflections 
6 6,;, and 6, for Input 1 


cf 


Vill 


eee @¢@¢ @ @ @ @ 


eee#e#eeee##ee#eee?@?#e#eee@e@.efhmUmMeMmCUmOrmUCUCOmUCUCOmUCPmUC mC OU OOhmhUCCrOhmhUCUc OhUlCUCrmhlUh 


Cr er 


e@eeeeeeeteeee#e e* # # #® #® @#® @# @® @# @#® @# @# # e@# ee ee # # # # @ @ 


YyH12.5 54 


Compensated X-29 Control Deflections 


ee ,, and ©. tor for Input 2 CTCEMC ss ss «ee eae 79 
Compensated X-29 Control Rates 
weer andsow mor EMNput 1 oo eee ct et ce ew eee can 
Compensated X=29 Control Rates 
or O¢, and 6, for Input 2 .......----2 eee ee eee eeee 83 


Limited-Performance X-29 H, Design 
Ea eree nS MEGS EEN eels sao 5 5 ss Stee wl ols SMe ee ee ete ee 86 


femeostmnecron |i, for y=2.62; 
Meni ea—PeTrLOrmanGe X=29 oes sos cess wee ec wee eee cece 87 


Sensitivity Function 8(s) for H, Solution, 
¥=2.62, Limited-Performance X-29 ...ceccvcccccsccces 88 


Complementary Sensitivity Function T(s) for H, 
Solution, y=2.62, Limited-Performance X-29 ......... 88 


Singular Value Plot I + G(s)F(s), 
mintriced-rerformance X—-29 .. cee ee www ew wee ec cet we tee 90 


Singular Value Plot I + (G(s)F(s))', 
mee ed — Pet fotemance X—29 ss. 6s sessile s s0 0 se ses ec sswes 90 


Singular Value Plot I + F(s)G(s), 


iermicea—=Performance X—-29 . oF. . ccc ccc ee cece e cc et tees 92 
Singular Value Plot I + (F(s)G(s))"', 
mmyced=Performance X—-29 ...cecccvccccccesscserevcecsses 92 


Limited-Performance X-29 a@ and gq Responses 


tem Nil (OCLC MNONSE OTN AEUC Te os ss es ss es ce ee ce ek eee ee ee 94 
Limited-Performance X-29 a and gq Responses 

re UG) sfc ee nels os «© 6 6 a sec no se ett tc ew te ce 95 
Limited-Performance X-29 Control Deflections 

or eecicleng PC MeNOUG To 6 ew es ess 5 se ct te ewe ee eee 96 
Limited-Performance X-29 Control Deflections 

meee  ONONO. LOR MAUL 2 as aie ee + switisis ee ee eee ness sie 98 
Limited-Performance X-29 Control Rates 

eo, emicleo., TOI OUT CI... see cs ewe ee eee 100 
Limited-Performance X~29 Control Rates 

me ON Oe Or EM pUE G2 8 se. we tw te ee OZ 


1x 


ACKNOWLEDGMENTS 


One justifiably experiences considerable satisfaction, as 
well as a tremendous relief, upon the completion of a thesis. 
Unfortunately, one is often too busy or too involved during 
the course of the thesis to thank those who contributed. But 
now as the pace slows and the final signatures are collected, 
it is time that I give the proper recognition and express my 
appreciation to the friends and associates who have provided 
advice, assistance, and direction. 

A special thank you goes to Professor Dan Collins, my 
thesis advisor, who suggested the topic of this thesis and 
provided the right amount of guidance and direction. His 
support was a key to its successful completion. 

A thank you is also extended to Professor Jim Hauser who 
served as the second reader to this document, and to Professor 
Art Schoenstadt of the NPS Mathematics Department. Professor 
Schoenstadt's generous help provided me an understanding of 
the mathematical theory underlying H, control synthesis. 

I also wish to express my gratitude to LCDR Tom McKannon 
who dealt with the "system" in procuring the computer hardware 
needed for this study, and to LT Russ Scott who unwittingly 
became a sounding board for my thoughts and ideas. Russ's 
only mistake was to occupy an adjacent computer work station 


throughout the course of this study. 


A most sincere thank you must go to June Melody, my 
friend, my love, my wife. Without her constant support, 
encouragement, and love this thesis could not have been 


realized. 


aa 











I. INTRODUCTION 


System performance, specifically a system's response to 
input commands, is of primary concern in the design of a 
feedback controller. Of equal importance in controller 
synthesis are stability margin and the remaining performance 
characterizations such as disturbance attenuation = and 
sensitivity. Stability margin is a relative measure of a 
control system's stability, i.e., defining the smallest 
perturbation or modeling error that will cause the control 
system to become unstable. Disturbance attenuation and 
sensitivity refer to a control system's ability to reduce or 
limit the effects of disturbances and plant variation, 
respectively, on the plant outputs. The analysis of single- 
input single-output (SISO) systems has the advantage of 
classical techniques, including Nyquist diagrams, Bode plots, 
and root locus plots, to measure the gain and phase margins 
and system performance. Attempts have been made to extend 
these classical SISO theories to the multi-input multi-output 
(MIMO) case. Often, however, these multivariable 
generalizations do not accurately represent the stability 
margin and performance characterizations of the MIMO system. 
A prime example is the multivariable generalization of the 


Nyquist Theorem. Although the stability of a MIMO system can 


be determined from a multivariable Nyquist diagram, the 
classical gain and phase margins are not meaningful. 

It has, therefore, become necessary to redefine the 
measures of multivariable stability margin and system 
performance. In recent years a number of control theorists, 
including Doyle, Safonov, MacFarlane, and Lehtomaki among 
others, have conducted considerable research on the analysis 
of multivariable system robustness [Refs. 1, 2, 3, 4, 5]. As 
a result of renewed concern with such multivariable feedback 
issues as stability margin, disturbance attenuation, and 
sensitivity, singular value Bode plots of return difference 
and loop gain matrices have emerged as useful frequency domain 
indicators of multivariable robustness [Refs. 1, 3}. The term 
robust or good robustness 1s used to describe a multivariable 
feedback control system characterized by a large stability 
margin, good disturbance attenuation, and low sensitivity 
[Ref. 3]. 

The development of H,, frequency-weighted linear quadratic 
guassian (LQG) (or H,), and LQG loop transfer recovery control 
theories, as well as the numerical optimization technique 
investigated by Gordon [Ref. 6], have made singular value loop 
shaping a reasonable if not routine procedure. By singular 
value loop shaping is meant the purposeful manipulation of a 
feedback control system's loop gains over ae specified 
frequency bandwidth in order to improve performance and 


stability margins. H, and H, theories provide direct, reliable 


techniques for synthesizing a controller which satisfies 
singular value loop shaping specifications. The LQG loop 
transfer recovery theory provides a less direct yet effective 
means of achieving singular value loop shaping requirements. 
The numerical optimization procedure is a relatively direct 
design method which manipulates the system feedback gains as 
design variables [Ref. 6:p. 4]. This manipulation is 
conducted such that singular value loop shaping requirements 
are met. [Ref. 7] 

Most published H, design examples have been textbook 
studies [Ref. 8]. It is the purpose of this thesis to assess 
the effectiveness of H, control theory in synthesizing a 
stabilizing controller for a reduced order, linearized model 
of the X-29 longitudinal dynamics. Chapter II discusses 
multivariable feedback properties, the role of the return 
difference matrix in evaluating these feedback properties, and 
Singular value loop shaping. H, control theory and design 
implementation are presented in Chapter III. The results of 
the H, synthesis of a controller for the X-29 longitudinal 
dynamics are discussed in Chapter IV. Conclusions are 


presented in the Chapter V. 


II. FEEDBACK PROPERTIES OF MULTIVARIABLE SYSTEMS 


The feedback properties of a linear system include 
stability and stability margin, sensitivity to plant and 
controller variations, and disturbance attenuation. These are 
the same system properties used in the Introduction to define 
the robustness of a multivariable feedback system, and can be 
altered only through the use of feedback. The 
characterization of a system's response to commands is omitted 
from this list of feedback properties since command response 
can be altered by prefiltering of the command signal. That 
is, feedback is not required for improving command response. 
This is not to imply that command response can not be 
controlled through feedback. However, purposeful manipulation 
of command response is best performed in the feedforward path. 
[tReft. 3] 

Throughout this thesis, vectors are denoted by bold lower 
case letters while matrices are indicated by bold upper case 
letters. 

A. RETURN DIFFERENCE MATRIX 

The return difference matrix is fundamental in the measure 
of a system's feedback properties and, therefore, to its 
robustness. Originally associated with SISO systems, the 
return difference concept has been extended to a MIMO 


generalization as discussed by Doyle and Safonov [Refs. 1, 3]. 


Consider the feedback configuration presented in Figure 
2.1, and let G(s) and F(s) be matrices of time invariant 
transfer functions for the system plant and controller, 


respectively. With the loop broken at the plant output, the 


transfer matrix 


I + G(s)F(s) (2.1) 





Figure 2.1 Return Difference Illustration 


is referred to as the output return difference matrix. This 
definition is more apparent if an external input vector u, is 
injected at point 1 of the loop. The transformed signal y 


returned at point 2 is 


Y(S) = -G(sS)F(s)u,(s), (2.2) 


and the difference between the injected input and the returned 


output vector is 
u4(s) = y(s) = [1% G(s) lise (2235 


If this loop is broken at the plant input, the resulting input 


return difference matrix is 
I + F(s)G(s). (2.4) 


The inverse-return difference matrices are defined for the 


output and input nodes as 


I + (G(s)F(s))"' (2.5) 
and 


Pee (Fr (Ss) G(s tie (2.6) 


respectively. Finally, the associated output and input loop 


gain matrices, also called return-ratio matrices, are 


G(S)F(s) (227) 
and 


F(S)G(S), (2.8) 
respectively. [Refs. 4, 9:p. 45] 


B. MATRIX NORMS AND SINGULAR VALUES 
The quantification of feedback qualities is necessary so 
that alternative feedback designs can be directly compared in 


the selection of an optimal design. To quantify the feedback 


properties of a system, a frequency-dependent, scalar measure 
of the return difference matrix size is required. For SISO 
systems the appropriate measure of the return difference 
matrix is its magnitude, i.e., |1 + g(jw)f(jw)]|. Eiisce is 
recognized as the distance to the critical point of the SISO 
Nyquist diagram for additive perturbation, and is a measure 
of the relative stability of the system. 

The notion of matrix size is extended to multivariable 
systems through matrix norms. The spectral norm, routinely 


used in the analysis of MIMO systems, is defined as 


|All, = max ,"(A"A) = o,,,(A) (2.9) 


max 
ot 


h 


where X, is the i eigenvalue of A’A, and A" refers to the 


complex conjugate transpose of matrix A. The singular values 
of a complex nxn matrix A, denoted o,(A), are the non-negative 
square roots of the eigenvalues of AA, that is 
_ Yo H 
a ee (AND) (2.10) 
Gmere 0, = 0, 2... 206, . The maximum and minimum singular 


n 


wearues of A are denoted by a... and o respectively. 


min? 


Some useful properties of matrix singular values are [Ref. 








S|: 
Ax 
ae o,., (A) = max 
x |x|]. 
| Ax || , 
e o...(A) = min 
x |x, 


3. oF, (A) SPX Gl =< oeey 
1 - 
4. if A’ exists, 0. () == — and 0,.,(A) = 
O max (A ) Onin(A ) 
5. o..,(@A) = |alo,.. (A) 


CncmmeO 2 (Amt Boe < «Gem, (Aa) ae (B) 
7. O,,,(AB) = 0 (Ajic= Gap 
8. max{o,.. (A), Ge,.(B)} = Geegiae)) 


< /2@max(o—. (A) , oNeGE oe 


Properties 1 and 2 are significant as they establish the 
Maximum and minimum gains of the matrix, respectively. The 
minimum singular value, property 2, is also used to measure 
a matrix's nearness to singularity. Finally, the 
multivariable generalization of return difference matrix size 


is 


o,{I + G(jw)F(j»)] (2.11) 
or 


o,{I + (G(jw)F(jw))'} (2.12) 


for additive and multiplicative perturbations, respectively. 
The significance of Equations (2.11) and (2.12) will be fully 


explored in the following section. 


C. MULTIVARIABLE ROBUSTNESS AND SINGULAR VALUE LOOP SHAPING 
Large loop gains, i1.e., G(sS)F(s) >> 1, are necessary to 


ensure good performance from a feedback control system. 


However, the need for ae system's tolerance to plant 
uncertainties restricts the frequency range over which the 
loop gains are permitted to be large. This tolerance to 
uncertainties, or modeling errors, iS a measure of a system's 
relative stability. Therefore, the design of a feedback 
control system requires a trade-off between performance and 
stability robustness. These concepts are fully developed in 
this section. 

Consider the multivariable feedback control system 
presented in Figure 2.2. Included are the transfer function 
matrices for the plant, G(s), and stabilizing controller, 
F(s), which are driven by the command, r, disturbance, d, and 
sensor noise, 7, vectors. in chrs seontiguration, d 4&5 
represented as the equivalent disturbance at the controlled 
or output variable vector, y [Ref. 9:pp. 168-174]. The 
feedback properties of this system, i.e., the multivariable 
stability margins and performance, can be measured using the 
closed loop transfer function matrices from the driving 
inputs, r, d, and 7n, to each of the outputs, y and the control 


vector, u. Namely, these matrices are 


S(s) = [I + G(s)F(s)]/ (2.13) 

T(s) = G(s)F(s){I + G(s)F(s)}"' =I - 8(s) (2% 
and 

N(s) = F(s)[I + G(s)F(s)]"' (2.35) 





Hy 
+ 
G 
MM <«<—np 
<i 


oe 


— aA 


>) F(S) >> G(s) > 


= <—7 


Figure 2.2 Feedback Control System 


where S(s), the sensitivity function, represents the transfer 
function matrix from d to y, T(s), the complememeaa 
sensitivity function, represents the transfer function matrix 
from r or -7 to y, and N(S) represents the transfer function 
matrix from r to u [Ref. 7]. The fundamental relationships 
between system performance and the above matrices can be 
realized by the following representations: 
1. > input-output behav tems) keri. 
¥(s) = G(s)F(s)[I + G(s)F(s)]7'(r-n) 
+ [I + G(stems ed (2.16) 


e= r-y 


[I + G(s)F(s)]'(r-d) + G(s)F(s)[I + G(s)F(s)]'7 


S(s)(r-d) + T(s)7 (2.17) 


10 


2. system sensitivity to variations [Ref. 9:p. 180]: 


(I + G(s)F(s) ]7'AG(s)F(s) [I + G(s)F(s)]/ 


AT (S) 


S8(s)AG(s)N(sS). (2.18) 
In Equation (2.18), AT(s) denotes the changes in the closed 
loop system as caused by plant variations or additive 


uncertainties, AG(s), where 
G(s) = AG(s) + G(s). (2.19) 


Equation (2.17) shows that system errors resulting from 
commands and disturbances can be made small by making the 
sensitivity function small or the return difference matrix, 
[I + G(s)F(s)], large. Equation (2.18) shows that closed loop 
changes, or loop sensitivity, is improved under the same 
conditions provided G(s) ~ G(s) [Ref. 9:p. 180}. Thus, the 
disturbance and sensitivity performance objectives of a 
multivariable feedback system can be represented in the 
following manner using the matrix singular values previously 


@iscussed: 


[I + G(jw)F(jw))7' < |w,' (jo) | (2.20) 


O max 


or equivalently, using singular value property 4, 


(I + G(jw)F(jw))] = |W, (je) | ZN 


O min 


where |W,(jw)| is a scalar function representing the desired 
disturbance attenuation and sensitivity factor. Equation 


(2.21) states that good feedback performance is achieved with 


1a 


large loop gains. This follows from the fact that the size 
of the return difference matrix approximates the loop gains 


whenever the loop gains are large: 


[I + G(s)F(s)] * (G(s) F(s) (2.22) 


Onin Onin 
fom 
(G(S)F(s)) >> 1. [Refs. 1, 7] 


Omin 


While large loop gains over a specified frequency range 
improve system performance with respect to disturbances and 
plant variations, they cause increased errors for large sensor 
noise inputs, 7. As seen in Equation (2.16), large loop gains 
result in G(s)F(s)[{I + G(s)F(s)]"', 1.e., T(s), near one. 
Therefore, sensor noise is passed through to the system 
outputs over the frequencies that the loop gains are large. 
It is obvious that a performance tradeoff must be made in the 
feedback design. [Refs. 1, 9:pp. 174-178] 

A more crucial limitation to high loop gains results from 
the need for tolerance to uncertainties. These uncertainties 
are due to deviations of the physical plant from its linear, 
time invariant model, and restrict the frequency range over 
which the loop gains may be large. The magnitude of the 
tolerance to these uncertainties is a measure of the stability 
margin enjoyed by the system. Representation of unstructured 
uncertainty, that uncertainty which is characterized solely 
by its magnitude, depends on the errors the model is expected 


to make. However, a common method of representation is the 
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multiplicative form. In this form, the true plant matrix 


G(jw) is given as 


G(jw) = G(jw) [I + E,(jo)] (2.23) 
with 

E,(jo) = G@'(jw)[G(jw) - G(jw))] 
and 


Cpe Teen) Wail) | Vw 2 0 


where E,(j#) is the relative (multiplicative) modeling error, 
G(jw) the nominal plant, and |W,(jw)| a scalar function which 
serves aS an upper bound on the relative error. That is, 
|W,(jw)| represents the size of the largest multiplicative 
modeling error anticipated for the respective frequency wo. 
The use of multiplicative error is preferable over absolute 
(additive) modeling error, denoted by Equation (2.19), since 
relative uncertainty applies equally to the compensated system 
G(s)F(s) and the raw plant G(s). This is not the case with 
absolute modeling error. [{Refs. 1, 5:pp. 73-85 and Refs. 7, 
meepp. 19-22) 

The properties of the bounding function |W,(jw)| are such 
that it is small (<< 1) at low frequencies and increases to 
one and above at higher frequencies (Figure 2.3). The maximum 
frequency W at which the bounding function is less than or 
equal to unity (the crossover frequency) is referred to as the 


multiplicative robust frequency, and is denoted by w As 


rm ° 


Will be discussed later, this frequency provides an upper 


3 


Magnitude 





Figure 2.3 Properties of the Relative Error Bounding Function 


constraint on the frequency band of the control system. A 
final comment concerning |W;(jw)| is that this bound assumes 
a worst case for the magnitude of the unstructured uncertainty 
that applies to all system loops. [Refs. 1, 7] 

Given the nominal plant model G(s), the performance 
objectives, and knowledge of the unstructured uncertainties, 
the problem is to synthesize a controller such that [Ref. 1} 


1. the feedback control system, G(s)F(s) [I + G(s)F(s)]" is 
stable, 


2. the stability margin is sufficient to cause stability 
of the actual feedback control system, G(S)F(S)3sm 
G(s) F(s))', for the largest anticipated relative error, 
ES) ecuing 


3. the performance objectives are achieved considering the 
langesteanticitpacred Ea (sie 
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These three conditions can be viewed as performance and 
stability margin specifications in terms of singular value 
loop shaping requirements [Ref. 7]. 

Condition 1 requires fulfillment of the multivariable 
generalization of the Nyquist Theorem [Ref. 5:p. 59 and Ref. 
6:p. 27]. The satisfaction of condition 2 is assured if and 


only if 


Onax( (JW) F (jw) [I + G(jw)F(jw))') < [Ws ' (30) | (2.24) 
or equivalently 


onfl + (G(jw)F(jw))'} > |W;(jo) |. (2.25) 


Equation (2.24) states that in order for a feedback system to 
be stability robust, its loop gains must be small when the 
unstructured uncertainty magnitudes are large, i.e., 

|Wz(jw)| >> 1. This follows from the fact that the 
complementary sensitivity function T(s) approximates the loop 


gain as the loop gain becomes small: 


c,.,{G(S)F(s) [I + G(s) F(s)]"'}) * o,,,.(G(s)F(s)) (2.26) 
mor 
(G(s)F(s)) << 1. 


O max 


(Refs. 1, 5:pp. 76-84 and Ref. 7] 


Equations (2.20) and (2.21) have already established the 
requirements for meeting the performance objectives as 
stipulated in condition 3. However, to account for the 


relative modeling error encountered at low frequencies, 


le 


|\W,(jw)| ais now considered to include the unstructured 
uncertainty for |W,(jw){ < 1. Therefore, restating the 


performance objective requirement: 


Omint(l + G(JW)F(jw)] 2 |W, (Je) | (2.27) 
for all © suche enae 
|W,(jw)| < 1 and o,,,(G(jw)F(jw)) >> 1 


min 


[Refs. 1, 7s 


The singular value loop shaping requirements established 
by the performance and stability bounds are depicted in Figure, 
ee The design of a MIMO feedback control system can be 


summarized as the use of high loop gains at low frequencies 


Ca (GF) 
O dB “ W 
PERFORMANCE STABILITY 
BOUND W, BOUND W,”' 
4 
7 
Onin (GF) é 


Figure 2.4 MIMO Singular Value Loop Shaping Requirements 
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where performance requirements are defined and the use of low 
loop gains at high frequencies where stability robustness 
constraints lie. While the low frequency bound on the loop 
gains is desirable to meet specified performance objectives, 
the upper frequency bound is a hard requirement. For if the 
uncertainties are unstructured and Equation (2.25) is 
violated, then there exists a modeling error E,(s) permitted 
by Equation (2.23) for which the system is unstable [{Ref. 1]. 
As alluded to earlier, the robust frequency w,, is an upper 
bound on the bandwidth w, of a multivariable control system 
(G(s)F(s)) >> 1. Thus 


where W, 1s the frequency range for One 


permitting the loop gain to be greater than unity above w,, 
violates the required condition for stability of the actual 
feedback system. The roll-off (attenuation) rate of the loop 


Meds, On. (G(S)F(s)) and o-. (G(s)F(s)), determine the severity 


min 
Geeene 0, constuaint. Large roll-off ratesspermit a wider w,. 
However, steep loop gain roll-offs are achieved at the expense 
baat (G(s) F(s))'] values when 


of small o._.. [I +G(s)F(s)]} and o 


min min 
the loop gains are approximately unity. This means both 
performance and stability margins are poor at the loop gain 
cross over frequency. This correlates with classical SISO 
analysis that suggests the slope of Bode magnitude plots be 
more gradual than -40 dB/decade to ensure an adequate 
eeiba tity margin [Ref. ll:p. 433]. [Refs. 1, 7] 


To better appreciate the design tradeoff required to 


achieve both performance objectives and stability robustness 


ds 


requirements, it is helpful to observe the algebraic relation- 


ship between o,. (I + G(s)F(s)] and o,. (I + (G(s) F(s))7'}. From 


min min 
Equations (2.25) and (2.27) it is seen that the bounds on low 
and high frequency unstructured uncertainty, as well as that 
on performance, become less restrictive if both of the above 
singular values can be made as large as_ possible. 
Unfortunately, these return difference quantities are related 


so both cannot be made large simultaneously: 


[I + G(s)F(s)] + (I +(G(s)F(s)) |} =H (2.28) 
equivalently 


(=) 4 (eee a. 


The bounded region of Figure 2.5 shows the allowable values 
of the minimum singular values of the return difference and 


the inverse return difference matrices. When the loop gains 











are large, I@., O-a(G(s)P(sjj) >> iat Lowecr 
Olt GCSt = ee or S(s) > 0 
while 
o.,(I + (G(s)F(s))7'] > 1 or MTs) > I 


Conversely, when the loop gains are small, 1.e., 0) (G(=) Riga 


< ] ation 





[I + (G(s)F(s))'] >> 1 or fs) > 0 


Onin 


while 








(I + G(s)F(s)] ea or S(s) >» I. 


Onin 


[Refs. 3, S52 ppaecl- 2.) 
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Onin 


[I + (GF)"'} 
A 


> 


[I + GF] 


Onin 
HeguLes2.o) Relacicnship Between o- [I + G(s)F(s)] and 
[I + (G(s)F(s))"'} 


Onin 


Finally, while o,. (I + G(s)F(s)] and o,, (I + (G(s) F(s))-'] 


are multivarible generalizations of SISO performance and 
stability margin concepts, it is important to note that these 
multivariable return difference singular values measure the 
uncertainty tolerances at the plant outputs only (see the 
development of Equations (2.1) and (2.4)). i Mee@iecer to 
measure the performance and relative stability at the plant 
{XI + F(s)G(s)]} and 


mputs, it is necessary to evaluate oa_.. 


[I + (F(s)G(s)°'}. 


Omin 


aS 


III. H, CONTROL DESIGN 


H, control theory provides a very powerful and direct, 
iterative procedure for the synthesis of a feedback controller 
designed to meet singular value loop shaping requirements. 
Its capabilities extend across the full range of performance 
and stability loop shaping objectives that can be formulated 
within the singular value Bode plot framework. The standard 
configuration of an H, problem is presented in Figure 3.1. 
The objective is to design a controller F(s) in state space 
form that stabilizes the augmented plant P(s) while satisfying 
specified performance and stability margin requirements as 


discussed in Chapter TPs kefs. & 2) 


P(s) 


F(s) < 


Figure 3.1 Standard H, Small Gain Problem 


This chapter briefly discusses the three major steps in 
the H, design procedure; the formulation of an H, control 
system problem, the selection of weightings to achieve design 
objectives, and the synthesis of an H, controller [Ref. 8]. 
In addition, a comparison is made among LQG, H,, and H, control 


syntheses. 


A. PROBLEM FORMULATION 

Consider the problem of designing a feedback controller 
F(s) for a multivariable plant with a transfer function matrix 
G(s) where G(s) = C(sI-a) 'B+D. The requirements for 
performance and stability robustness are defined by Equations 
(2.27) and (2.25), respectively, in terms of the system's 
return difference matrices. These requirements are redefined 
here in terms of the sensitivity S(s) and complementary 


sensitivity T(s) functions as 


1 
— 2 17, (34) | W < W (oll 
of eG) a}2)) 

and 
PCs | wee (30) | W > U, Gr) 


where w, and wo, are the respective O @B crossover frequencies 
op yW, and W,"' and y is a positive scalar constant. Gamma (7) 
is added to facilitate iterative attenuation of the 


sensitivity function during design implementation. [{Ref. 12] 
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Within the» framework sof H, controle) seneoGma eed. 
requirements of Equations (3.1) and (3.2) can be combined into 


a single infinity norm specification as 


Ve eten= 2 (3.3) 


where 


(3.4) 


eis) is the closed loop transfer function of the augmented 
feedback control system shown in Figure 3.1 from input) Ujee 
output y,. The H, norm of a transfer function matrix K(S)ie 
denoted in terms of its frequency-dependent singular values 
o.(jW) as 

IKI, = 7 Snax (K(j¥)) (3.5) 


where "sup" or "supreme" refers to the least upper bound [Ref. 


7]. Therefore, the H, small gain problem |7,,,,l|, < 1 states 
that the maximum singular value of Triui(S), as definediaam 
Equation (2.9), is to be less than or equal to one for all 
frequencies w. As a consequence of singular value property 
Syl, 

a A A 

O max Ss Max{On,,(A), Crane s O max . 
2 B B 


Ze 


Equation (3.4) approximates the requirements of Equations 
(3.1) and (3.2) to within a factor of /2, or 3 dB [Refs. 7, 
a2 | . 

The development of Tyiuy Can be shown with the plant G(s) 
augmented by the weighting functions YW,, W,, and W, as 
presented in Figure 3.2. As evident from the figure, YW, 
weights the error signal e, W, weights the control inputs u, 
and W, weights the plant outputs y,.. The area enclosed by 
asterisks is P(s), the augmented plant, previously shown in 


Figure 3.1 where 





Pi, | Pi 
P(s) = a (S50! 
Po, | Poo 
or 
yw, | -71W,¢ 
o| WwW, 
P(s) = (3.7) 
0 | WG 


as 


Y) P,, 








; (3.8) 
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Figure 3.2 Augmented Plant 


P,,U, + PU, 


P,,u, + P,U. (3.9) 


u, may be expressed as 


Fy> (3.10) 


be written as 


P,,u, + P,Fy>. (3.11) 
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Putting y, in terms of u,, F, and P;; and substituting into 
Equation (3.11) gives 
Y 


“4 
Dee =—eP. + P.F(I — P,.F) P,,- (3.12) 
u 





Substitution of the appropriate P,;,; elements from Equation 


(3.7) into Equation (3.12) results in 


1W,S 
Tut = WN |. (3.13) 


WT 


Ignoring the W, weighting function gives Equation (3.4). The 
fmmeeenction can be used to weight or penalize the control 


input u to the plant G(s). 


The augmented plant P(s) has a state space realization of 


[Ref. 7] 


PS eee oD. Dd, |. (3.14) 


By putting each weighting function and the plant G(s) from 


Figure 3.2 into standard state space form, i.e., 
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Ax + Bu 


“ 
ll 


Cx + Du, 


he! 
I 


Equation (3.14) may be represented as 











A. J 8 o 0 B, 
~ByiS, Ay, 90 0 Bi), B,,,D, 
0 oO aA, oO 0 Bs 
ye) = B30, 0 0 Aus 0 By3D, 
“Dy Cy, 0 0 Dy ~D,,D, 
0 0. Omer 0 0 De 
Df, 0 0 Cy 0 DD, 
-C, 0 0 0 I -D, 


(3.15) 


As is apparent from Equation (3.4), the closed loop 
transfer function Duis) 1S partitioned into submatrices 
representing the performance and stability robustness 
requirements. The weighting functions yW, and W, serve as low- 
pass or high-pass filters in order to emphasize or de- 
emphasize, aS appropriate, the maximum singular values of 
Tyiui(s). These weighting functions are applied over the range 
of frequencies necessary to achieve the desired performance 
objectives while meeting design constraints [Ref. 8]. 

Therefore, the control problem defined by the singular 


value loop shaping requirements of Equations (3.1) and (3.2) 


Ae 


is realized with the solution of the H, small gain problem 
(Figure 3.1) [Ref. 7]: 
Given a transfer function matrix P(s), find a stabilizing 
controller F(s) such that the closed loop transfer 


muUNGtION Tyiu1 (S) is internally stable and its infinity 
norm is less than or equal to one: 


IT yuu (50) I, < 1. (3.3) 
The compensated system shown in Figure 3.1 is said to be 
internally stable if the A matrix of the augmented plant P(s) 
is stable. In other words, with the external input u, equal 
to zero, all states of the P(s) and F(s) transfer function 
matrices will asymptotically go to zero for any initial 
condition [Ref. 8]. As indicated by Francis [Ref. 13:pp.34- 
35], a sufficient condition to stabilize P(s) is for F(s) to 
stabilize P,,(s). 

In general, the solution to the H, small gain problem is 
not unique since a number of stabilizing controllers will 
satisfy Equation (3.3). Conversely, a solution to the small 
gain problem does not exist for every P(s) due to violation 
of well-posedness conditions or infeasible singular value Bode 
plot specifications. In order for an H, stabilizing 
Sempensator to exist, it is necessary for (A,B,) and (A,B,) 
to be stabilizable and (C,,A) and (C,,A) to be detectable. 
Additionally, for the infinity norm specification of Equation 
(3.3) to be realized, Tyiui(S) must have no eigenvalues on the 


imaginary jw axis. If such poles exist then |T,,,l], > 1. 


Finally, it is necessary for the D,, and D,,' submatrices of 
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Equation (3.14) to have full column rank. The Di> submatuae 
physically weights the plant's control inputs. That is, the 
matrix [Ds Dis is comparable to the LQG control weighting 
matrix R of the full state cost function. A full column rank 
D,. matrix ensures all control inputs are weighted, thereby 
eliminating infinite impulses. The W, weighting function can 
be used to ensure D,, 1s full column rank as is evident in the 
augmented plant P(s) state space realization, i.e., Equation 
(3.15). [Refs. 7a 

The software program used to solve the H, control problem 
outlined in the following chapter is hinf which is included 
in the Matlab Robust-Control Toolbox package. Hinf employs 
a variant of a two Riccati formula of Doyle et al. The 
advantage of using hinf to solve the small gain problem is 
that the two Riccati formula eliminates the lengthy and 


numerically sensitive model reduction work characteristic of 


earlier algorithms. The H, controller produced by hinf has 
the same number of states as the augmented plant P(s). [Ref. 
wv) 


B. SELECTION OF WEIGHTS 

The weighting functions yW, and W, are transfer function 
matrices whose diagonal elements are frequency-dependent 
constants. The size of these weighting matrices is consistent 
with the number of plant output states. As the weighting 


functions act as filters, their design parameters include 


2° 


gain, corner frequency, attenuation, and order. The goal in 
the selection of these design parameters is to achieve small 
singular values of the sensitivity function S(s), i.e., large 
loop gains, over the broadest possible bandwidth subject to 
the constraints imposed by stability robustness, i.e., the W, 
design constraint. [Refs. 8, 12] 

The weighting functions are directly associated with the 
physical characteristics of the feedback control system being 
designed. These characteristics include time domain 
specifications (rise time, settling time, and overshoot), 
performance bounds, and command response in the case of YW,, 
Hegmrcudewor thevcontrol inputs in the case of W,, and largest 
Pmetcipacted relative modeling error E,(s) in the case of W,. 
It should be understood that the resultant H, controller is 
Only optimal with respect to the yW,, W,, and W, weights 
selected. It is necessary, therefore, that the designer have 
a good insight into the physical capabilities of the system, 
including stability robustness constraints, to ensure that a 
reasonable H, problem is posed. Postlethwaite [Ref. 8] 
suggests selecting a fixed weighting W, representing the 
stability robustness constraints, and as large a weighting 
yW, aS possible representing system performance. The 
objective is to make the largest singular value of S(s) as 
small as possible over some operating frequency band by making 
yW, as large as possible over the same frequency band. This 


is done while satisfying Equation (3.3). The following 


Zo 


algorithm is recommended in choosing appropriate weights for 
the H, control problem {Refs. 7, 8]: 


1. Select the weight W, to represent the stability 
robustness (design) constraints. 


2. Select the weight W, and set y=1 to represent an inasges 
guess of the performance objectives over the desired 
operating bandwidth. The yW, 0 dB crossover frequency 
™. must be less than the W, crossover frequency \@ayas 
order that the stability robustness constraint, Equation 
(3.2), is not violated. 


3. Solve the small gain problem in search of a stabilizing 
Controller that satisfies Equation (3.3), 2-em 


Ital. <2. 
4. If |T,,ll. < 1, select @ larger 7 and réepeam ete rues 


5. Continue increasing y until a stabilizing controller no 
longer exists that satisfies Equation (3.3), 1.ée., 
NT yau1 tle > 1 or the stability robustness constramam 
Equation (3.2), is violated. 


6. If it is necessary to make the D,, submatrix full column 
rank, choose a W, weighting with an invertible D matrix, 
1.@., Wo(s)=el where € 185 any non=Zzeroenuneem A W, 


weighting can also be used to regulate the magnitude of 
the control energy input to the system plant G(s). As 
the magnitude of € is increased relative to 7, a larger 
weighting is placed upon the control inputs u. This 
results ina smaller control energy input to the plant 
G(s). Decreasing the magnitude of € relative to y has 
the opposite effect. 
Cc. CONTROL SYNTHESIS 
Once the design requirements are specified in terms of the 
weighting functions yW,, W,, and W,;, the H, controllermam 
calculated within a computer aided design environment. The 


procedure involves an iterative process accomplished by 


increasing the coefficient y. Increasing y, or decreasing 
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Y , results in a smaller DC gain and higher wo, 


for (7 - 
Therefore, the iterative increase of y results in a continual 
suppression of the sensitivity function S(s) and a widening 
ef the control bandwidth 4,. This iteration is continued 


until the sensitivity function S(s) is forced against its 


upper constraint of (yW,) 7! Bom UemGmeine) 3.3). Due to Che 
N 7 
PERFORMANCE 4 Yo. 
BOUND (YW,) yi aN 
O dB : = W 
Lo STABILITY 
, BOUND W,_ 
we 
O max (S) ue 
a 
O (T) rb 


- and W; 


Figure 3.3 Singular Value Bode Plot of S, T, W, 
relationship between the sensitivity function S(s) and 
complementary sensitivity function T(s), 1.e., Equation 
meeeeo), T(s) goes to I as S(s) goes to QO. Therefore, 
increasing y forces the complementary sensitivity function 


T(S) against its upper constraint of Ws £ Ole W>W. ChaGur,emsn:3 ) 


oA 


[ReEs a2). Figure 3.4 presents a flow chart of the H, 
iterative procedure [Ref. 7]. Augmentation of the plant G(s) 
with the weighting functions and presentation of the necessary 
Singular value Bode plots are performed by separate M-files 


within the Matlab Robust-Control Toolbox package. 


D. COMPARISON OF LQG, H,, AND H, CONTROL SYNTHESES 

The intent of the following discussion is to provide an 
appreciation for the advantages of using H, control synthesis 
over the use of other modern methodologies such as LQG and H, 
norm (also known as frequency-weighted LQG). This end will 
be accomplished by revealing the disadvantages and short- 
comings of these latter two synthesis methods. Only a cursory 
discussion of the LQG and H, methods will be presented. Those 
unfamiliar with these methods are referred to the listed 
references for a more thorough explanation. 

In the following discussion, reference to linear quadratic 
(LQ) regulators refers to those regulators with full state 
feedback. LQG refers to regulators that have output feedback 
and use a Kalman filter to provide state estimates for 
feedback. LQG problems involve the solution of two Riccati 
equations, one to solve for the optimal state feedback gain 
and a second to solve for the optimal filter gain. 

It is well documented [Ref. 5:pp. 191-211] that linear 
quadratic (LQ) regulators possess guaranteed minimum gain and 


phase margins of -6 dB to © and +60 degrees, respectively. 
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This is true for a diagonal control weighting matrix R. These 
guaranteed stability margins imply "desirable" loop gain 
properties which, as discussed in Chapter II, define system 
performance and stability robustness. Unfortunately, these 
LQ stability margins and associated loop gain properties do 
not necessarily apply at the plant inputs and outputs of LOG 
regulators. As discussed by Lehtomaki [Ref. 5:pp. 217-226], 


the Kalman filter is the dual of the LQ regulator and does 


possess the above guaranteed stability margins. The Kalman 
filter, therefore, possesses the desirable loop- gain 
properties. 


The standard LQG control system block diagram is presented 
in Figure 325. The points marked 2 and 3 represent the 
feedback control system input F(S)G(s) and output G(s)F(s) 
loop gain matrices, respectively. It is at these locations 
that the guaranteed stability margins of the LQ problem are 
desired. However, it is at the loop gain matrices of points 
1 and 4, inside the Kalman filter loop, that the guaranteed 
stability robustness properties apply. This limitation to LQG 
robustness is caused by the Kalman filter's inability to 
account for modeling errors. That is, the Kalman filter 
models the nominal plant and not the actual plant as given in 
Equation (2.23). [{Refs. 1, Ssppee2zi7—-27c8 

There are, however, two dual procedures that will 
adequately recover the loop gain properties at the plant input 


(point 2) or output (point 3) to approximate the loop gain 
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functions at points 1 and 4, respectively. Each of these two 
procedures, referred to as "full-state loop transfer 
recovery," permits the recovery of desirable performance and 
stability robustness properties at the plant input or output 
aS appropriate. A limitation to these recovery procedures is 
that they work only for a minimum phase plant G(s). An 
additional limitation is the loop gain recovery applies only 
at the input or output of the physical plant butG ee 
necessarily at both. As performance and stability robustness 
qualities are important at both the plant input and output, 
some design tradeoffs must be made. Finally, the regulator 
or Kalman filter gains necessary to effect the loop gain 
recovery at the plant input and output, respectively, are 
often very large and physically impractical. [Refs. 1, 5:pp. 
226-229 and Ref. 7:pp. RR42-RR44} 

A primary difference between LQG and H, control syntheses 
is the absence of frequency-dependent shaping of the loop 
gains with the LQG methodology. Thus, H, control synthesis 
provides the designer more control over the shaping of 
performance and stability robustness attributes due to the 
inclusion of the weighting functions yW, and W,. H, control 
synthesis is a frequency-weighted LQG process that, like H,, 
uses the frequency-dependent  yW,, W,, and W, weighting 
functions and augmented plant P(s). In the case of Chemie 


control problem, however, the objective is to calculate a 


a6 


stabilizing controller F(s) that minimizes the H, norm of the 
closed loop transfer function T,,,,(s), that is {Ref. 7] 


min|T ll, 
where 


- en 
oe ll, = | a ste (Tau (j®)) ys aw }. (3.16) 
Ban as 


Using the Matlab h2lqg M-file, the H, norm control problem 
is solved as a conventional LQG problem involving two Riccati 
equation solutions, one each for the full state regulator and 
Kalman filter optimal gains. The H, control problem 
formulation, weight selection, and synthesis format is 
adentical to the H, control problem. H, synthesis is often 
used with H, synthesis as a first cut to determine the levels 
of performance achievable. Figure 3.6 presents a flow chart 
Seeene combined H,/H, synthesis process. 

An interesting feature of the H, two Riccati equation 
solution is that the resultant controller displays a 
mearation structure similar to that observed in LQG or H, 
problems. That is, the algebraic Riccati equations provide 
solutions to what can be considered a suboptimal state 
feedback gain and a suboptimal filter gain. [{Ref. 7] 

Although H, synthesis shares the singular value loop 
Shaping capabilities of the H, procedure, the H, designed 
controller usually does not match the robustness’ levels 


attained by the H, design. That is, an H, designed controller 
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will, in general, provide a larger bandwidth of operation and 
greater disturbance and uncertainty attenuation within the 


stipulated stability constraints. 


oe 


IV. X-29 H, CONTROLLER SYNTHESIS 


A challenging application of H, control theory is the 
synthesis of a stabilizing controller for the longitudinal 
dynamics of the X-29. The X-29 is a technology demonstrator 
with a unique forward swept wing design that offers the 
advantages of improved maneuverability, better low speed 
handling, and reduced stall speeds. The X-29's longitudinal 
dynamics are designed with 35 per cent negative static 
stability margin [Ref. 15]. Unlike the advanced fighter model 
with an unstable phugoid mode used by Safonov [Ref. 7] to 
demonstrate H, and H, methodology, the X-29 has an unstable 
short period mode, i.e., a real pole on the positive axis. 

The X-29 controller synthesis was performed using Pro- 
Matlab and the Matlab Robust-Control Tool Box software. These 
application packages were run on a Sun 386i work station. The 
script files specifically written or modified for this problem 
are listed in Appendix A. 

This chapter discusses the X-29 state space model, the H, 
design objectives and specifications, the controller 
synthesis, and the design results including the aircraft's 


longitudinal responses to test inputs. 


40 


A. X-29 MODEL DESCRIPTION 

The X-29 longitudinal dynamics model is that of the 
aircraft's analog reversion mode with the aircraft trimmed at 
-5 mach, 30,000 feet. An 83rd order model was reduced to a 
14 state model that includes a short period approximation of 
the aircraft longitudinal dynamics, vertical velocity w and 
pitch rate q, and fourth order actuator dynamics for each of 
the three longitudinal control surfaces, i.e., the canards, 
flaps, and strakes. Eliminated from the 83rd order model were 
the flexible mode dynamics, aerodynamic lag terms, sensor 
dynamics, and notch filter. 

Figure 4.1 presents the physical configuration of the open 
loop actuator/aircraft dynamics model. Some of the actuator 
gains shown in Figure 4.1 may have changed in the current 
aircraft configuration. For the purposes of this study, two 
separate commands, r, and rz, are input to the three control 
eueface actuators with r, controlling the canards and r, 
controlling the flaps and strakes. Although not truly 
representative of the X-29, this configuration provides 
multiple, independently controlled surfaces representative of 
a super-maneuverable aircraft. As will be seen, this 
configuration results in the synthesis of advanced control 
modes which are characteristic of super-maneuverable aircraft 
designs. The control inputs to the aircraft dynamics are the 
canards 6,, the flaps 6,, the strakes 6, and their respective 


first and second derivatives. The measured outputs are the 
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two aircraft states, wandq. Thus, the model has two inputs, 
two outputs, and 14 states. 

The 14th order model was scaled in order to improve the 
numerical conditioning of its state space representation. In 
the scaling process the w state was transformed to angle-of- 
attack a, i.e., @ = w/u, where u, is the initial forward 
velocity, and the units of the actuator third derivative 
states were transformed from rad/sec to 1e+04 rad/sec’. This 
scaling was effective in reducing the condition number of the 
system's A matrix from an order of magnitude of Ok ae OmelOns 

The state space realization of the resultant 14 state 
linear model G(s) = C(sI - A) 'B+D is presented in Appendix B. 
The order of the state variables, their description, and 
respective units are listed in Table 4.1. Finally, the open 
loop poles of the X-29 model are listed in Table 4.2. Note 


that the unstable short period mode has a real pole at 1.9550. 


B. DESIGN OBJECTIVES 

The singular value plot of the xX-29 plant G(s) is 
presented in Figure 4.2 where the solid curve is g,.,(G(jw)) 
and the dashed curve is o,.,(G(jw)). The uncompensated X-29 
model possesses poor disturbance attenuation, high sensitivity 
to variations and modeling errors, and a small control 
bandwidth w, as evidenced in Figure 4.2 by the small loop 
gains at the lower frequencies. These performance 


characterizations are to be improved by suppressing the 
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TABLE 4.1 


ORDERED LIST OF THE UNCOMPENSATED X=29 MODEEDSS A 


State 


Description 


angle-of-attack 
pitch rate 

Canard Control inpue 
flap control Uampme 
strake control input 
CanaragecontEoteaace 
flap controMecace 
strake control rate 
canard cConmtrel accel 
flap Control accer. 
strake control accel. 
canard control jerk 
flap control jerk 


strake control jerk 
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Units 


rad 
rad/sec 
rad 
rad 
rad 
rad/sec 
rad/sec 


rad/sec 


rad/sec® 


rad/sec* 


rad/sec* 


1e+04 rad/sec° 


1e+04 rad/sec” 


le+04 rad/sec° 


TABLE 4.2 


UNCOMPENSATED X-29 OPEN LOOP POLES 


2.3201le+021 


ete 


-2.2746e+02 
-1.4491e+02 
-1.4455e+02 


1.9550e+00 
-1.0031le+02 
-2.7155et+00 
-5.2506e+01 
-5.2518e+01 
-5.0067et+01 
-2.0172e+0O1 
-2.0115e+01 


4.8410e+011 
+ 4.8255e+0O11 


a 


X-29 SV PLOT (G) 
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Singular Value Plot of the X-29 Plant G(s) 


Figure 4.2 
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sensitivity function singular values o,(8(jW)) as much as 
possible, i.e., making the loop gains as large as possible, 
over as wide a bandwidth as possible. These performance 
objectives must be performed within the bounds of the system's 
Stability constraints. In this problem, it is necessary to 
attenuate the closed loop singular values of the complementary 
sensitivity function o,(T(jw#)) by 20 dB at frequencies beyond 
100 rad/sec, and to exhibit a second order roll-off beyond 100 
rad/sec. These stability requirements ensure that the X-29 
system has sufficient stability margin to tolerate modeling 
errors or loop transfer function variations which could arise 
from the unmodeled flexible modes. These flexible modes are 
observed in the Bode plot of the 83rd order X-29 model (Figure 
4.3) at the frequencies 100 rad/sec to 250 rad/sec. The 
second order roll-off also closely matches that of the open 
loop plant (Figure 4.2). 

The following (YW,)'(s), W.(s), and W; '(s) weighting 
functions were selected to meet the above performance 
objectives and stability constraints: 


iz .01(100s + 1) 








(yW,) '(s) = y * — ——_ ne 
-O1ls + (Zax 2) 

A B 

W5(S) = = -, 001s 
Cc D (4 xX 4) 

2 1000 

W, (S) = 5 * I ; (4.1) 

Ss (2. Xas2) 
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Figure 4.3 Open Loop Bode Plot of X-29 83rd Order Model 


A plot of the H, W, '(s) and W, '(s) weighting functions is shown 
MpFigure 4.4. The W.(s) weighting is included to ensure the 
D,, submatrix of the augmented plant P(s) has full column rank 
equation (3.15)). This weighting function penalizes the 
control u input to the X-29 plant G(s) as shown in Figure 3.2. 
As the W.(s) weighting function is not proper, it has no state 
Space realization. However, the term W,(s)G(s) seen in the 
matrix representation of the augmented plant P(s) (Equation 
(3.7)) is proper and can be realized in the required state 
space form. The augx29.m script file listed in Appendix A 
performs the matrix calculations necessary for this 


realization. The results of this matrix manipulation are 
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X-29 Design Specifications 
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Figure 4.4 xX-29 H, Design Specifications 


placed into the D,,C, and D,D, terms of the augmented plant's 
C, and D,, submatrices, respectively, by the augx29pl.m script 
file (Appendix A). 

The resultant X-29 augmented plant P(s) is an 18th order 
system with W,(s) and W,(s) each adding two states to the 
X-29 plant G(s). The W;(s) function does not add states to 
G(s) as a state space realization of this function does not 
exist. Using the two Riccati solution, H, and H, synthesized 
controllers are the same size as the augmented plant P(s). 


As such, an 18th order X-29 controller is expected. 
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CONTROLLER SYNTHESIS 


Ce. 


The H,/H, synthesis process outlined by Figure 3.6 was 


followed with the H, solution used as an initial indication of 


achievable performance levels. 


To ensure a well-posed H, 


(s) weighting 


a 


1 


the upper corner frequency of the W 


problen, 


denominator a 


function was removed by making the W, '(s) 


This ensured that the D,, submatrix of the 


constant value. 


augmented plant P(s) is 0 as required by H, control theory. 


A plot of the H, design specifications is presented in Figure 
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Using the x29h2.m and x29hinf.m script files listed in 
Appendix A, solutions to the H, and the H, small gain problems 
were obtained for increasing values of y until: 


1. For the H, solution, the cost function |T,,,,ll, reached 
the "all pass limit", i.e., O GB. 


2. For the H, solution, no stabilizing controller satisfied 
the H, small gain problem for a larger value of jy, 1.e., 
no solution existed for a larger j. 
Figures 4.6 and 4.7 are singular value plots of the maximum 
(solid curve) and minimum (dashed curve) singular values of 
the H, cost function |T,,,,ll, for y values of one and 6.7, 
respectively. Figure 4.8 is an identical plot for the H, cost 
function |T,,,,l, for a y of 12.5. Figures 4.6 and 4.7 Show 
that as y increases from one to a maximum of 6.7 in the H, 
solution, the maximum singular value of |T,,,,|, increases to 
0 dB. The H, solution pushes both ||T,,,,l, singular values to 
within 2 @B of the "all pass limit" (Figure 4.7). However, 
with a maximum y of 12.5, the H, solution pushes the |[T,,,,|l. 
Singular values to within .5 dB of the "all pass limit" 
(Figure 4.8). The significantly larger value of y in the H, 
solution indicates that a higher level of performance is 
achieved with the H, solution compared to that reached by the 
H, solution. 
Figures 4.9 through 4.11 are singular value plots of the 
sensitivity function 8(s) and (yW,) '(s) weighting function for 
the H, solution with y values of one and 6.7, and Chesae 


solution with a y of 12.5, respectively. The dotted curves 
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Figure 4.11 Sensitivity Function 8(s) for H, Solution, y=12.5 


are o,..(8(jW)) while the dashed curves are o,,.(S(jW)). 
Figures 4.12 through 4.14 are singular value plots of the 
complementary sensitivity function T(s) and W, '(s) weighting 
function for the H, solution with 7 values Of one and 69 7)ame 
H, solution with a y of 12.5, respectively. Here, the dashed 
curves are o,,,(T(jw)) while the dotted curves are o,, (T(jW)). 
As y¥ is increased from one to 12.5, the sensitivity function 
S(s) is incrementally suppressed by the (yW,) '(s) weighting 
function, and the complementary sensitivity function T(s) is 
pushed toward the stability constraint, 1.e., the S(s) and 


T(s) singular values are forced against their respective upper 


lalalne (yW,) '(s) and W; '(s). Comparing Figures 4.10 and 4.11, 
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it is seen that o,(8(jwW)) is suppressed more by the H, solution 
meaning the H, compensated X-29 is characterized by larger 
disturbance attenuation, lower sensitivity to plant variations 
and modeling errors, and a wider control bandwidth 4a,. 
Inspection of Figures 4.13 and 4.14 shows that the closed loop 
bandwidth of the H, solution is larger, reaching nearly 30 
rad/sec. This indicates the H, compensated X-29 is a more 
responsive aircraft than the H, compensated aircraft. 

As anticipated, the H, solution to the small gain problem 
results in an 18th order controller. To "clean up" the 
solution, a minimal realization was performed using the Matlab 
minreal M-file to eliminate two uncontrollable and 


unobservable states. The minreal M-file finds the similarity 
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transformations of a system's state space realization such 
that the A, B, and C matrices are put into staircase 
controllability and observability forms. Minreal then removes 
the isolated uncontrollable or unobservable states from the 
system model. Next, the minimal 16th order controller was 
balanced using the Matlab obalreal M-file to improve the 
numerical conditioning of the controller model. The obalreal 
M-file finds a similarity transformation such that the 
mappings from inputs to states and states to outputs are 
balanced while preserving the input-output relationships. 
Finally, a fast transient pole at -1.6262e+05 was removed 
using the Matlab modred M-file. Modred eliminates specified 
states from a state space model while preserving the system's 
input-output relationships. Each of the 16 minimal states was 
individually eliminated until the state that removed the fast 
transient pole was identified. The resultant 15th order 
minimal controller was used in plotting the H, cost function 
ia lle: sensitivity function S(s), and complementary 
sensitivity function T(s) results presented in Figures 4.8, 
4.11, and 4.14, respectively. 

However, further reduction of the H, controller size is 
desirable. Therefore, attempts were made to reduce the 15th 
order controller using the Schur additive error model 
reduction method [Ref. 7]. This model reduction method allows 
the size of the reduced order model to be selected. Additive 


modeling error E, (s) is defined as 
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E,(S) = F(s) - F(s) (4.2) 


where F(s) and F(s) are the true transfer function matrix and 
its reduced model, respectively. To ensure the additive 
modeling error E, will not destabilize the closed loop systen, 
the control bandwidth w, of the reduced model must be less 


than the additive robust frequency @.., 1-e-, 


where 
W., = max{w | o,, (F(jw)) 2 0,,,(E,(jw)) }- (4.3) 
This condition for system stability assumes 


o (ERG) < oe Ceaiens for wv < ee (4.4) 


ra 


That is, given E, is the only information available about the 


modeling error, the additive robust frequency w,, 1s an upper 


a 
bound on the bandwidth of a multivariable control system 
without violating the required condition for stability, i.e., 
Equation (2.25). [Ref. 7] 

Figure 4.15 compares the singular value plot of the 15th 
order controller with that of a selected eighth order 
controller. The upper (solid) curve 1s a plot of ot jae 
and Onax (F (jw) ) for the 15th order and 8th order controllers, 
respectively, while the lower (dashed) curves are plots of 
(F(jw)) and o,, (F(jw) as indicated. Figure 4.16 displays 


Omin 


58 


ORDER X-29 CONTROLLERS 


Bth & 15th 





FREQUENCY - rad/sec 


Figure 4.15 Comparison of 8th Order and 10th Order 


29 Controllers 
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8th ORDER X-29 CONTROLLER ERROR & ERROR BOUND 
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Figure 4.16 8th Order Controller Additive Error BE, 


a9 


the singular values of the eighth order controller's additive 
error where the dashed and dotted curves are ga... (E,(jW)) and 
oi,(B,(j¥)), respectively. The solid curve is an upper bound 
on g,.,(E,(j¥)). As shown in Figure 4.15, the additive Tope 


frequency W.., 1.e., where 


Onax(E(J0)) = O,;,(F(j)), 


1s approximately 2 rad/sec. The control bandwidth wo, of over 
200 rad/sec, defined at the 3 dB point, is clearly greater 
than the additive robust frequency @... As the closed iiga 
feedback control system can be destabilized by the additive 
modeling error §E,, the eighth order controller 15 (noe 
suitable alternative for the 15th order controller. 

However, examination of Figure 4.17, which compares the 
Singular value plots of the 15th order and a selected 10th 
order controller, shows the two controllers to be well 
matched. Figure 4.18 shows the oa... (E,(j¥)) and a. (E,(jw) ) 
plots of the 10th order controller's additive error, i.e., the 
dashed and dotted curves, respectively. Again, the solid 
curve 1s an upper bound on @a,.,(E,(jw)). As can be seen from 
Figures 4.17 and 4.18, oa. (E,(j¥)}) does not interseen 
Oa eGI@))s) within the frequency range of the graph, §i-c yj 


> en rad/sec. Clearly, the control bandwidth ®,; of Che D@mg 


order controller, shown just below 200 rad/sec, is less than 
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Figure 4.17 Comparison of 10th Order and 15th Order 


29 Controllers 
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Figure 4.18 10th Ordér Controller Additive Error E, 
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the additive robust frequency w and the 10th order 


ra‘ 
controller is a suitable alternative to the 25th opeaew 


cOntise: ber. 


D. DESIGN RESULTS OF THE H, CONTROLLER 

The closed loop configuration of the H, compensated X-29 
is presented in Figure 4.19 where F(s) and G(s) are the 15th 
order controller and 14th order X-29 plant transfer function 
matrices, respectively. The output vector y is made up of the 
alrcraft states a and q, and the command vector r is composed 
of separate command elements r, and r,. Unlike the open wags 
actuator/aircraft dynamics model presented in Figure 4.1, fF, 
and r, are reference commands to the controlled outputs, @ 
and q. This is a result of the controller being placed in 


series with the xX-29 plant in the feedforward loop and the 





Figure 4.19 Feedback Configuration for H, Compensated X-29 
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negative unity feedback. Thus, the closed loop, compensated 
X-29 model has two inputs, two outputs, and 29 states. As 
will be seen, this closed loop configuration still provides 
the multiple, independently controlled surfaces observed with 
the open loop, uncompensated X-29 plant G(s). 

The state space realization of the 29th order, closed loop 
model is presented in Appendix B. A balanced realization was 
performed separately on the controller F(s) and the X-29 plant 
G(s) to improve their numerical conditioning. As a result of 
this balancing, the internal structures of F(s) and G(s) were 
altered making identification of the individual states 
difficult. The poles of the closed loop model are listed in 
Table 4.3. It is interesting to note that the unstable short 
period pole of the open loop system iS mirrored into the left- 
half plane of the closed loop system, i.e., -1.9550. This is 
not a coincidence as an identical occurrence is observed in 
Safonov's advanced fighter example presented in Ref. 7. In 
Safonov's example, the unstable phugoid poles of the advanced 
fighter's open loop model are mirrored into the left-half 
plane of the compensated, closed loop model. This) Mirror 
imaging can represent a basic limitation to the system's 
performance if, as in the case of the compensated X-29, this 
is the dominate pole. 

As discussed in Chapter II, the singular values of the 
return difference and inverse-return difference matrices 


quantify a system's feedback properties. In the following 
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TABLE 


H, COMPENSATED X-29 


4.3 


CLOSED LOOP POLES 


-4.1327e+02 
-2.2745e+02 + 2.3201e+021 
-1.3068e+02 + 4.6111e+011i 
-5.3305e+01 + 8.9700e+011 
-1.4491e+02 
-1.4452e+02 
-~1.3014e+02 
-1.3877e+01 + 5.9243e+011 
-9.9794e+01 
-5.2545e+01 + 4.8359e+01i 
-5.2503e+01 + 4.8301le+01i 
-7.4132e+01 
-1.9550e+00 
-2./155e+00 
-2.0379e+01 + 2.1564e+01i1 
-2.0578e+01 + 1.8907e+01i 
-4.9199e+01 + 6.5316e+001 
-4.2465e+01 
-2.0184e+01 
-2.0110e+01 


paragraphs the feedback properties of the H, compensated X-29 
will be measured uSing singular value plots of its return 
difference matrices. 

The singular value plots of the uncompensated and 


compensated X-29 output, return difference matrices are 


presented in Figures 4.20 and 4.21, respectively. In these 


figures, as with all the return difference matrix plots, the 


upper curves are 0... while the lower curves are 0... Recaam 


X 


from Chapter II that the output return difference matrix 
I + G(s)F(s) (also referred to as the output, additive return 


difference matrix) is the inverse of the sensitivity function 
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S(s), and that the minimum singular value of this return 
difference matrix approximates the loop gains whenever the 
loop gains are large, i.e., Equation (2.22). Given this, 
Figure 4.20 indicates the uncompensated X-29 possesses small 
loop gains along with the corresponding traits of poor 
Gisturbance attenuation and high sensitivity to plant 
variations and modeling errors. Figure 4.21 shows how the H, 
synthesized controller has markedly improved the X-29 
performance properties. The large loop gains indicate good 
disturbance attenuation and low sensitivity to uncertainties 
over a control bandwidth 0, of approximately 9 rad/sec. 
However, the dip of the singular values below the O dB line 
indicates that performance is lacking near the 0O dB crossover 
frequency. This 1S probably caused by the steep roll-off 
(-40 dB/decade) designed into the complementary sensitivity 
function. Finally, the 9 rad/sec control bandwidth w, of the 
compensated X-29 is less than the multiplicative robust 
frequency w,, of 30 rad/sec ensuring that the required 
condition for stability (Equation (2.25)) 1s not viovatedmae 
the, X-29 plantwoutpuc. 

The singular value plot of the output, inverse-return 


difference matrix I + (G(s) F(s))' 


(also referred to as the 
output, multiplicative return difference matrix) for the 
compensated X-29 is shown in Figure 4.22. The minimum 


singular value of this matrixjei-e.) sae ba (G(jo) (F(jo)) 
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Figure 4.22 Singular Value Plot PeGisiens).) ee 
H, Compensated X-29 


can provide a measure of the aircraft's gain and phase margins 
with respect to multiplicative modeling errors using the 
universal gain and phase margin curve (Figure 4.23) [Ref. 
Sepp. 54-55]. As shown in the universal curve, a minimum 
Singular value of one, i.e., O dB, provides gain and phase 
margins of -6 dB to infinity and +60 deg, respectively, or 
that stability margin guaranteed by the LQ regulator problem. 
It is seen from Figure 4.23 that any singular value less than 
one is associated oe suboptimal gain and phase margins. 
Pemurning Co Figure 4.22, it is seen that o, (I + 
Maley) F(jv))~'} drops to approximately -2 dB at frequencies 
between 1 rad/sec and 20 rad/sec. Entering the vertical axis 
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Figure 4.23 Universal Gain and Phase Margin Curve 


of the universal curve at a singular value of .79, i.e., -2 
aB, the gain and phase margins of the compensated X-29 near 
the 0 dB crossover frequency are -14 a@B to +5 aB and 147edeqe 
respectively. This is more stable than the -8 @B to +4 QB, 
+35 deg gain and phase margins typically designed into a 
fighter aircraft. 

The singular value plots of the input additive and input 
multiplicative return difference matrices are shown in Figures 
4.24 and 4.25, respectively. These singular value plots show 
that the H, solution to the X-29 small gain problem does not 
satisfy the performance objectives or the stability 


constraints at the X-29 plant input. An exceptionally poor 
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FREQUENCY - rad/sec 
Singular Value Plot I+(F(s)G(s))"', 


H, Compensated X-29 


Figure 4.25 


disturbance attenuation and low sensitivity to uncertainties 


are indicated by a small a,, [I + F(}#)G(j0)], especian 
between the frequencies of 2 rad/sec and 300 rad/sec (Figure 
4.24) 2) A Oe aie (F(jw)G(jw))7'] of -11 dB equates to a gain 
and phase margin of -3 dB to 2 dB and +15 deg, respectively, 
over much of the bandwidth between 2 rad/sec and 60 rad/sec 


(Figure 4.25). Further, since 
Onin(I + (F(jW)G(jw))'] < [Ws (jo) |, 


the required condition for stability, Equation (2.25), is 
violated meaning a relative modeling error exists that can 
destabilize the aircraft at its inputs. It can be concluded 
that H, control synthesis does not guarantee that the 
performance and stability robustness requirements of Equations 
(2.27) and (2.25) will be satisfied at the plant G(s) inputs. 

A .01745 rad (1 deg) step input was applied separately for 
one second to each of the two reference commands, and the 
aircraft responses and control deflections plotted. Following 
standard convention, a positive a or q corresponds to a nose- 
up response, and a positive control deflection is trailing 
edge down for all three control surfaces. Positive canard 
deflection 6,, negative flap deflection 6,, and negative strake 
deflection 6, each induce a nose-up response, i.e., positive 
OO and *¢- 

The closed loop, compensated X-29 model (Figure 4.19) 


exhibits precision flight path control modes as a result of 
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the multiple, independently controlled surface configuration. 
Figure 4.26 presents a graphic representation of these 
precision control modes in terms of angle-of-attack a, pitch 
attitude @, flight path angle y, and the aircraft principal 


and stability axes x, and x,, respectively. (Note that the 


p 
flight path angle y is in no way related to the y constant 
used in the H, and H, design specifications. ) The three 
precision longitudinal modes observed are [Ref. 3}: 

1. Vertical Translation: The aircraft vertical velocity 
is controlled at a constant @ by varying a, i.e., the 
aircraft flight path angle y or velocity vector is 
controlled while x, remains fixed. 

2. Direct Lift Control: The aircraft flight path angle y 
is controlled at a constant a by varying @, i.e., the 
aircraft flight path angle y or velocity vector remains 
along the aircraft stability axis x, as x, rotates. 

fe Pitch Pointing: The aircraft pitch attitude @ is 
controlled at a constant flight path angle jy, i.e., the 
aircraft flight path angle y or velocity vector remains 
fixed while x, rotates (@ = a). 

The q@ and q responses of the compensated X-29 are 
presented for inputs 1 and 2 in Figures 4.27 and 4.28, 
respectively. The compensated X-29 responds to input 1 (r,) 
With a positive a and negligible change in gq, i.e., order of 
magnitude is 10°, (Figure 4.27) which is the vertical 
translation mode discussed above. That 1s, input one 
decouples q and @ from a. The compensated aircraft exhibits 
a fast response to input 1 with an a rise time of .125 sec. 
The direct lift control mode is effected by input 2 (r,) as 
shown in Figure 4.28. The aircraft responds to input 2 with 


a positive q while its a response is negligible, i.e., order 
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Figure 4.28 Compensated X-29 a and q Responses to Input 2 


a> 


of magnitude is 10°. As awalen input 1, a decoupling of eae 
aircraft responses, a and gq, is observed with input 2. The 
aircraft q response is equally fast with a rise time of 
approximately .095 sec. Figures 4.29 and 4.30 show the 
compensated X-29 control surface deflections 6,, 6,, and 6, for 
inputs 1 and 2, respectively. From these figures it is 
determined that the flap and strake control inputs, 6, "angie 
respectively, control the aircraft's a response while the 
canard 6, controls the q (and @) response. A simultaneous 
injection of inputs 1 and 2 is necessary to effect the pitch 
pointing precision control mode. 

The control deflections 6., 6,, and 6, for the) Vere 
translation mode (input 1) have peak magnitudes between 1.4 
rad and 3 rad which exceed the xX-29's control surface 
deflection limits of; 

1. canards (leading edge): 30 deg up / 60 deg down, 

2. flaps (trailing edge): 10 deg up / 25 deg down, 

3. strakes (trailing edge): 30 deg up and down. 
The control deflections are significantly less for the direct 
lift control mode (input 2) with peak magnitudes between .15 
rad and .6 rad, and are more closely aligned with the control 
surface deflection limits listed above. A similar difference 
in control rates between inputs 1 and 2 is observed in Figures 
4.31 and 4.32 which present the compensated X-29 control rates 
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) 6,, and 6, for inputs 1 and 2, respectively. The peas 
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Figure 4.29 Compensated X-29 Control Deflections 
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Figure 4.29 (Cont) Compensated X-29 Control Deflections 
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Figure 4.30 (Cont) Compensated X-29 Control Deflections 
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Figure 4.31 (Cont) Compensated X-29 Control Rates 
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Figure 4.32 


83 


X-29 DSDOT FOR 0.01745 rad / 1 sec STEP (INPUT 2) 


30 


20 


— 
© 


radians/second 
°o 





TIME — sec 


Figure 4.32 (Cont) Compensated X-29 Control Rates 
6 6,, and 6, for Input 2 


cl 
control rates for both inputs 1 and 2, ranging from 
approximately 8 rad/sec to over 100 rad/sec, greatly exceed 
the X-29 actuator, minimum design requirements of; 

1. canards: 1.75 rad/sec (100 deg/sec), 

2e4 filaps: .87 rad/sec (50 deg/sec), 


3. strakes: .52 rad/sec (30 deg/sec). 


E. LIMITED-PERFORMANCE H, CONTROLLER DESIGN 
The H, synthesis of a stabilizing controller for the X-29 
was reworked to bring the control surface deflections and 


control rates more into line with physical capabilities. A 
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secondary objective was to retain the system robustness 
achieved with the initial H, controller. To accomplish these 
objectives, a greater weighting or penalty was applied to the 
control input vector u of Figure 3.2. This was effected by 
increasing the magnitude of the W,(s) ¢€ term to .025. 
Additionally, the upper corner frequency of the W, '(s) 
weighting function was moved from 100 rad/sec to 10 rad/sec 
to restrain the sensitivity function overshoot near the O dB 
crossover frequency. The weighting function assignments for 


the limited-performance controller design are: 


Peel 200s + 1) 





(yW,) '(s) = y! * —————_ «I 
is sell Cae) 
A B 
NS). eee = -.025 * I 
Cc D (4 x 4) 
Z 1000 
ves) — *1 (4.5) 
S (Quay 


A plot of the W, '(s) and W, '(s) weighting functions is shown 
in Figure 4.33. 

To avoid confusion, the initial H, controller design will 
be referred to as the optimal-performance case. This is not 
meant to imply that the solution is the result of a cost 
function minimization as with the LQ regulator problem. 


Rather it implies the compensated X-29 performance is optimal 
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Figure 4.33 Limited-Performance X-29 H, Design 
Specifications 
with respect to the design specifications. The reworked 


controller design will be termed the limited-performance case. 

With a maximum achievable y of 2.62, the limited- 
performance H, solution is only able to push the Eales cost 
function singular values to within 5 a@aB of the "all pass 
limit" (Figure 4.34) as compared to .5 @aB for the optimal- 
performance H, solution (Figure 4.8). The singular value 
plots of the sensitivity function 8(s) and the complementary 


sensitivity function T(s) for the limited-performance X-29 
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Figure 4.34 H, Cost Function |T,,,,l., for y=2.62, 
Limited—Performance X~-29 


are presented in Figures 4.35 and 4.36, respectively. These 
figures are from the H, solution for a 7 of 2.62. Comparison 
of Figures 4.35 and 4.36 with Figures 4.11 and 4.14 shows that 
the feedback properties of the limited-performance X-29 are 
not as robust as those of the optimal-performance aircraft. 
The limited-performance X-29 is characterized by smaller 
disturbance attenuation, larger sensitivity to plant 
variations and modeling errors, a smaller control bandwidth, 
and a smaller closed loop bandwidth. 

The state space realization of the 29th order, closed loop 
limited-performance model is presented in Appendix B, and the 
Closed loop poles are listed in Table 4.4. Note that, once 
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TABLE 4.4 


CLOSED LOOP POLES OF THE LIMITED-PERFORMANCE X-29 


—-2.2747e+02 + 2.3201le+02i 
-1.3612e+02 
~-1.4415e+02 
-1.4494e+02 
~1.4476e+02 
-1.0023e+02 
-5.1688e+01 
-7.2313e+t0Ol 
-~3.7606e+01 
~5.1716e+01 
-5.3298e+01 
~5.2512e+01 
-5.0505e+01 
-3.6342e+01 
-1.1889e+01 + 1.2160e+01i1 
-1.9539e+00 

-2.7204e+00 

-3.7372e+00 

-9.6752e+00 

-2.0387e+01 + 1.1180e+001 
-2.1180e+01 


I+ 


7.7445e+01i 


5.2777e+011 
5.0560e+01i 
4.7220e+011i 
4.8306e+011 


Pecitesl eit 


again, the unstable short period pole of the open loop system 
was mirrored into the left-half plane of the closed loop 
system, l.e., -1.9539. 

Review of the output return difference matrices confirms 
the above characterizations of the limited-performance 
aircraft (Figures 4.37 and 4.38). The singular value plot of 
O,,,(L + G(jJW)F(jW)] (Figure 4.37) clearly shows a smaller loop 
gain as compared to that of the optimal-performance X-29 (Fig- 
ure 4.21). Referring to Figure 4.38, o,, [I + (G(jw)F(jw))'} 
drops to approximately -3.45 dB, a .67 singular value, between 
1 rad/sec and 10 rad/sec. Using the universal gain and phase 


margin curve of Figure 4.23, this equates to a gain and phase 
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FREQUENCY - rad/sec 


Singular Value Plot I + (G(s) F(s))"', 


Limited-Performance X-29 


Figure 4.38 
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margin of -10 daB to +5 dB and +439 deg, respectively. 
Therefore, the stability margins of the limited-performance 
X-29, with respect to multiplicative modeling errors, is 
closer to the values desired of a fighter aircraft, i,e., -8 
dB to +4 dB, +35 deg gain and phase margin, respectively. 
The singular value plot of the input, additive return 
difference matrix (Figure 4.39) shows that the limited- 
performance H, solution also did not satisfy the desired 
performance objectives at the X-29 plant inputs. That is, 
G.,LL + F(j#)G(j#)] is small, especially between .8 rad/sec 
and 30 rad/sec. However, the singular value plot of the 
input, inverse-return difference matrix (Figure 4.40) shows 
that the required stability constraints are met at the X-29 


plant inputs, l.e., 


Onintl + (F(jw)G(jw))') > |Wz(je) |. 


Although the stability constraints are satisfied, the relative 
stability with respect to multiplicative modeling errors is 
very low between .7 rad/sec and 15 rad/sec. This is evidenced 
Pyeene large sub-0 dB values of o,,[I + (F(jw)G(jw))7'). A 
minimum singular value of -13 dB at 4 rad/sec equates to a 
gain and phase margin of -2 dB to 2 dB and +12 deg, 
respectively. 

Identical step inputs were applied So the limited- 


performance xX-29 model as were previously applied to the 
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X-29 SV PLOT (I + FG) 
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FREQUENCY - rad/sec 


Singular Value Plot I + F(s)G(s), 


Limited-Performance X-29 


Figure 4.39 


X-29 SV PLOT (I + inv(FG)) 
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Limited-Performance X 


=29 


Figure 4.40 
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optimal-performance model. Figures 4.41 and 4.42 show the a 
and q responses of the limited-performance X-29 for inputs 1 
and 2, respectively. As with the optimal-performance X-29, 
the limited-performance X-29 responds to input 1 (r,) with a 
positive a while input 2 (r,) effects a positive q response. 
However, the decoupling of the a and gq responses is not as 
pronounced as observed in the optimal-performance X-29. This 
indicates that the precision flight path modes are not 
achieved in the limited-performance X-29 to the extent 
observed in the optimal-performance case. The limited- 
performance X-29 responses to the step inputs are slower with 
rise times of .5 sec and .22 sec for a@ and q, respectively. 
This is a result of the smaller closed loop bandwidth observed 
in the complementary sensitivity function. 

The control deflections 6,, 6,, and 6, for the limited- 
performance X-29 (Figures 4.43 and 4.44) are considerably 
smaller than observed for the optimal-performance case. The 
peak magnitudes of the deflections for inputs 1 and 2 vary 
from .1 rad to .33 rad and, with the exception of negative 
flap deflection for input 1, are well within the X-29 control 
surface deflection limits. It is also seen that the input 1 
and 2 deflections are more closely matched. Similar 
observations are made for the control rates Sn Oe and 6. for 


inputs 1 and 2 shown in Figures 4.45 and 4.46, respectively. 
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X-29 alpha RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 1) 


0.02 


0.01 


radians 


0.005 





TIME - sec 


X-29 q RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 1) 
0.015 


0.01 


0.005 


radians/second 


-0.005 


-0.01 





TIME - sec 


Figure 4.41 Limited-Performance X-29 a and q Responses 
to Inpuiees 
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x10-3 X-29 alpha RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 2) 


radians 





TIME — sec 


X-29 q RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 2) 
0.025 


0.02 


0.015 


0.01 


radians/second 


0.005 





2 


TIME —- sec 


Figure 4.42 Limited-Performance X-29 a and q Responses 
to Input 2 
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X-29 DC FOR 0.01745 rad / 1 sec STEP (INPUT 1) 
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i 

TIME — sec 
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Figure 4.43 Limited-Performance X-29 Control Deflections 
6 6,, and 6, for Input 1 


ac? 
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X-29 DS FOR 0.01745 rad / 1 sec STEP (INPUT 1) 


radians 





TIME - sec 


Figure 4.43 (Cont) Limited-Performance X-29 Control Deflections 
6 6,, and 6, for Input 1 


cc? 


or) 


X-29 DC FOR 0.01745 rad / 1 sec STEP (INPUT 2) 


0.08 


0.06 


radians 


0.04 


2 


TIME -— seco 


X-29 DF FOR 0.01745 rad / 1 sec STEP (INPUT 2) 


0.25 
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radians 


0.05 





0 0.5 1 1.5 2 250 3 3.5 4 


TIME -— sec 


Figure 4.44 Limited-Performance X-29 Control Deflections 
6 6,;, and 6. for Inpwem 


c? 
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X-29 DS FOR 0.01745 rad / 1 sec STEP (INPUT 2) 


0.25 


radians 


0.05}. 


0 : : 2 





TIME - sec 


Figure 4.44 (Cont) Limited-Performance X-29 Control 
Deflections 6, 6,, and 6. for Input 2 
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X-29 DCDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1) 
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TIME — sec 


X-29 DFDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1) 
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Figure 4.45 lLimited-Performance X-29 Control Rates 


6 6,, and b. for Input 1 
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X-29 DSDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1) 


radians/second 
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Figure 4.45 (Cont) Limited-Performance X-29 Control Rates 
52 6,, and 6, for Input 1 
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X-29 DCDOT FOR 0.01745 rad / 1 sec STEP (INPUT 2) 


0.6 


0.4 


0.2 


radians/second 
° 


=—(re 





TIME - sec 


X-29 DFDOT FOR 0.01745 rad / 1 sec STEP (INPUT 2) 


radians/second 
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TIME — sec 


Figure 4.46 Limited-Performance X-29 Control Rates 


bo, Se and 5. for iInpiee 


OZ 


X-29 DSDOT FOR 0.01745 rad / 1 sec STEP (INPUT 2) 


radians/second 





TIME - sec 


Figure 4.46 (Cont) Limited-Performance X-29 Control 
Rates Soaec,, and 62 for Input 2 


That is, the peak control rates of the limited-performance 
X-29 range from .4 rad/sec to approximately 7 rad/sec which, 
with the exception of 6. for input 2, exceed the actuators' 
capabilities. 

Although the primary objective of the limited-performance 
H, controller was achieved, i1.e., the control deflections and 
control rates were reduced, it is obvious from the results 
that the performance robustness of the compensated X-29 was 
reduced. Considering the results of the optimal-performance 
and limited-performance H, solutions, there is a direct 


relationship between the level of system performance attained 
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and the weighting of the control inputs u (Figure 3.2). Small 
control weightings permit a greater influx of energy into the 
plant G(s) and improve system performance, i.e., larger loop 
gains and a larger control bandwidth. The constraint for this 
particular example is the physical plant's inability to 
accommodate the energy levels necessary to attain the desired 
levels of performance. That is, actuator performance is 
inadequate for the compensated X-29 to achieve the performance 
levels specified in the H, small gain problem. The limited- 
performance solution works around this constraint by 
penalizing the control input u more heavily and accepting a 


reduced system performance. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


H, control theory provided a systematic and effective 
procedure to synthesize a stabilizing controller for the 
statically unstable X-29. The framework of H, control theory 
ensured that the performance and stability design 
specifications, i.e., the  yW, and W, weightings, were 
satisfied. As a result, the H, optimal-performance X-29 
demonstrates (at the plant outputs) good disturbance 
attenuation, low sensitivity to plant variations and modeling 
errors over a control bandwidth of 9 rad/sec, and sufficient 
stability margin to tolerate perturbations from the unmodeled 
flexible mode dynamics. 

Additionally, H, control theory demonstrated an ability to 
cope with a true multivariable design problem such as the 
multiple, independently controlled surfaces characteristic of 
a super-maneuverable aircraft. The H, optimal-performance 
X-29 demonstrates a quick response to command inputs along 
with a decoupling of its angle-of-attack and pitch modes. 
This latter attribute is necessary to effect the advanced 
control modes in which aircraft flight path and pitch attitude 
are independently controlled. However, the results also show 
that the control surface actuators are saturated by the high 
input energy necessary to realize the performance 


specifications. That is, the control surface deflections and 
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control rates required to effect the desired performance 
attributes exceed the present actuators' performance 
capabilities. 

It was shown that a work-around to the problem of actuator 
saturation is an increased penalty on the control input u. 
The magnitude of the control surface deflections and control 
rates for the limited-performance X-29 model are more in line 
with physical capabilities at the expense of reduced 
performance robustness. Therefore, it 1S apparent that 
physical plant capabilities are an additional constraint in 
achieving the desired levels of performance, and must be 
included in the H, problem formulation. 

While H, control theory ensured satisfaction of the design 
specifications at the X-29 plant outputs, the theory did not 
ensure these specifications were met at the plant inputs. 
This occurred since the feedback property specifications 
are formulated in terms of the sensitivity functienm 
(I + G(s)F(s))7' and complementary sensitivity function 
G(s)F(s) (I + G(s)F(s))°' which are structured in terms of the 
plant outputs. 

It is recommended that the synthesis of an H, stabilizing 
controller for the X-29 be repeated with actuator performance 
considerations included in the H, problem formulation. This 
can be accomplished by including the control inputs 677g 
and 6esand/ors@he coneroleraees sa o ana 8. in the output 


vector y, (Figure 3.1). The objective of this formulation @e 
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to design a stabilizing controller that will not saturate the 
control surface actuators, yet have the compensated X-29 match 
the performance and stability robustness achieved with the 
optimal-performance case. 

Additionally, it 1s recommended that a second order W,(s) 
weighting function be considered in the H, problem formulation 
in order to further suppress the sensitivity function S(s) 


at frequencies immediately below its 0 dB crossover frequency. 


Oe] 


%*x29h2.m 
diary x29h2.dat 
format short e 


APPENDIX A 


X-29 MATLAB SCRIPT FILES 


disp(' a) 
disp(' This script file is designed to solve the H2 control ') 
disp(' problem for the X-29. The 14th order FDLTI model, in state ') 
disp(' space form, is that of the X-29 aircraft and actuator dynamics.') 
disp(' Two states are those of the aircraft dynamics, i.e., alpha ') 
disp(' and q. The remaining 12 states are the dynamics of the ') 
disp(' three, fourth order actuators, i.e., the canard, flaperon, and ') 
disp(' strake actuators. The order of the unbalanced states is as ') 
disp(' follows:') 
disp(' alpha, gq, dc, df, ds, dcdot, dfdot, dsdot;edcdeiaorcnas 
disp(' dfdbldot, dsdbldot, dcetrpldot, dftrpldet, astrplaae® 
disp Given the open loop transfer function G(s)=Cinv(IS-A)B+t+D, a ') 
disp(' Stabilizing controller F(s) will be found such that the H2 norm') 
disp(' of Tylul is minimized. ‘') 
disp(' H2 optimal control synthesis is performed to determine attainable’) 
disp(' performance levels. Once completed, an Hinf control ') 
disp(' Synthesis is performed.') 
disp(' ") 
pause 
cle 
disp(' The scaled X-29 aircraft and actuator state space representation’ ) 
disp(' a 
ag=[—.4181d+00 .9960d+00 -.2269d-01 -.1213d+00 -.1948d-01 -.9493d=025 
-4427d-04 -.6712d-04 .1451d-05 -.2162d-04 -.3540d-05 0.0 OJ OMNGaa 
~-5474da+01 -—.3424d+00 -2585d+01 -.1386d+01 -.1058d+01 .3898d=-025e 
-~.1164dad-01 -.6397d-02 -.2509d-03 -.5362d-03 -—.2912d-03550) 03) 
0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 90.0 0.0 0:0) (O00 3 o0o 
0.0 0.0 0.0 0.0 0.0 90.0 1.0 0.0 0.0 0.0 0°30) 50 (GG 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0050 GOs 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0 0.0 9020 30)? 
0.0 0.0 0.0 0.0 0.0 0.0 9.0 0.0 0.0 1:50 (O20 550° G0] 
0.0 0.0 90.0 0.0 0.0 0.0 0.0 0.0 0.0 O [0°31 0356 oe ee 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.061045 0°03 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 90.0 0.0 0.0 0.0 0.0 1.06+O2RaeS 
0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0O°0 a 
0.0 0.0 -.1479d+04 0.0 0.0 ~-1143d+03 0.0 0.0 -.2529d+01 Ostia 
0.0 -.2697d+03 0.0 0.0; 
0.0 0.0 0.0 -.1491d+04 0.0 0.0 -.1149d+03 0.0 0.0 — 325364) 
0.02 °0.0) = 270 lass O..0.: 
0.0 0.0 0.0 0.0 =-.5302d+05 0.0 0.0 -— J 18l6a7i4 7 
-.1790d+02 0°50 0.0 =. 0053d7037 
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bg=(0.0 0.0; 
O20 0.0; 
©.0 O60; 
e.0 O20; 
O70 0.0; 
e270 0.0; 
oO 0.0; 
e.0 0.0; 
oO 0.0: 
m20 0.0; 
c.0 640.0; 
-1479d+04 0.0; 
0.0 .1491d+04; 
0.0 .5308d+05] 
aw—it.0 0.0 0.0 90.0 90.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 
moe l.0 0.0 0.0 0.0 60.,0mNO.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0} 
dg=zeros (2) 
pause 
disp(' ") 
disp(' Balanced realization of the X-29 state space representation' ) 
disp(' a) 
{agbl,bgbl,cgbl,g,t]jJ=obalreal(ag,bg,cgq) 
ag=agbl; bg=bgbl; cg=cgbl; 
disp(' ~) 
disp(' Calculate the poles and transmission zeros of the open loop plant') 
disp(' ") 
poleg=eig(ag), tzerog=tzero(ag,bg,cg,dg) 
disp(' Ya 
disp(' e 
disp(' Determine the condition number of ag') 
disp(' ") 
disp(' ) 
condag=cond(ag), rcondag=rcond(ag) 
disp(' i) 
disp(' i) 
pause 
disp(' ) 
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disp(' << Design Specifications >> ') 

disp(' ry 

disp(' 1). Robustness Spec. : -40 dB roll-off, -20 db @ 100 Rad/Sec.") 
disp(' Associated Weighting: ') 

disp(' ") 

disp(' aaa 1000 ') 

disp(' W3(8) =  ------ * I (fixd)') 

disp(' 2) 

disp(' 8 2x2') 

disp(' “) 

disp(' ) 

disp(' 2). Performance Spec.: minimizing the sensitivity function’) 
disp(' as much as possible.') 

disp(' Associated Weighting: ') 

disp(' " 

disp(' '*) 

disp(' ae -1 (100s + 1) i) 

disp(' W1(s) = Gam * errrercre---- ~*~ io) 

disp(' 100 2x2) 

disp(' a) 

disp(' where "Gam" in this design is iteratively updated from 1') 
w=logspace(-3,3,100); 

k=1000; mn=[(2 2); tau=0.0; 

nuw3i = [0.0 k]); dnw3i = [1.0 O Oj; 

Svw3i = bode(nuw3i,dnw3i,w); svw3i = 20*1l0g10(svw3i); 

nuwli = [100.0 1.0]; dnwli =[{0 100.0]; 

Svwli = bode(nuwli,dnwli,w); svwli = 20*1log10(svwli); 

aw2=-.00l*eye(2); bw2=zeros(2); cw2=zeros(2); dw2=-.001l*eye(2); 


disp(' ") 

disp(' -) 

disp(' (strike a key to see the plot of the weightings ...)') 
pause 

semilogx(w, Svwli,w,Svw3i) 

Grid 


title('X-29 Design Specifications') 
xlabel('Frequency - Rad/Sec' ) 
ylabel('1/W1l & 1/W3 - dB') 
text(.01,0,'1/W1(s)') 

text (.5, 100,.° Wwsacs)\ 

meta x29specs 

pause 

ele 
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disp(' 
disp(' ) 
disp(' 
disp(' =) 
disp(' 
disp(' ") 
disp(' 
disp(' i 
disp(' 
disp(' 
disp(' ") 
disp(' 
disp(' 
disp(' iy) 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' i) 
disp(' ") 
disp(' ) 
disp(' 
pause 

cic 

disp(' ") 
disp(' 2) 
disp(' 
disp(' iy) 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' ) 
disp(' a) 
disp(' 
pause 

elc 

disp(' -) 
disp(' ) 
disp(' ") 


where 


a) 


zx kek kK ke ke kK Ke Ke Ke KK Ke Ke KK Ke kK Ke KK Ke Ke Ke Ke kK kK kK kK 


[Step 1}. Do plant augmentation (run AUGMENT.M or 


+ € € % & € HH F BF ¥ 


<< Problem Formulation >>') 


ioe 


ne 


Mirek 


F(s) 


i 
J 
Tylul = | 
| 
I 


[Step 2}. 


[Step 3]. 


{Step 4]. 


i i a a a a a a a a a a 


Gam*W1*(I + GF) |} = 


Lee. ') 


Gam*Wl penalizing error signal "e"') 


W3 penalizing plant output "y"') 


< 1,") 


Zita) 


=5) 7) 
Gam * Wl * S 
w3 * (I - S) 


! 
1 
“1} 


W3*GF*(I + GF) |') 


(strike a key to continue .. 


<< DESIGN PROCEDURE >>') 


AUGX29.M) 


new "Gam" 


Do H2 synthesis (run H2LQG.M) 


Redo the plant augmentation for a 


and rerun H2LQG.M 


Redo the plant augmentation for a 


higher "Gam" then run HINF.M 


(strike a key to continue ... 


aL 


and find a stabilizing controller F(s) such that the H2 norm') 
of TF Tylul is minimized, 


1") 
I 
I 


‘) 


=) 0) 


)Fee) 


Form an augmented plant P(s) with these two weighting functions:') 


Ee) 
er) 
-) 
a) 
a) 
E. ) 
» 
a) 
») 
- 
e 


x") 


disp(' 
disp(' 
disp(' 
pack 


Assign the cost coefficients "Gam" with Gam=1 ') 
. 


serving as the baseline design .. 


saa) 


gama=(1.0,6.7]; 
ngama=length(gama) ; 
for i=l:ngama 
Gam=gama(1,i) 


k=1000 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
augx29 
disp(' 
disp(' 
disp(' 
pause 
cle 
disp(' 
disp (' 
$disp(' 
disp(' 
disp(' 


") 


> 
") 


(strike a key to continue ...)') 


2 
") 
Do state space balancing on the augmented plant if needed') 
") 
") 


{abal,bbal,cbal,g,t])=obalreal(A,([Bl B2)],(C1;C2)) 


$A=abal, 


disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 


Bl=bbal(:,1:2), 
2) 

") 

The transmission zeros, poles and condition number of the augmented’ ) 
plant follow. In addition, determine if (A,Bl) & (A,B2) are ') 
stabilizable and if (C1,A) & (C2,A) are detectable.') 

3 

) 


B2=bbal(:,3:4), Cl=cbal(1:6,:), C2=cbhal(7:6)7) 


tzeroaug=tzero(A, (Bl B2),(C1;C2],(D11 D12;D21 D22)) 
poleaugA=eig(A) 

condaugA=cond (A) 

rcondaugA=rcond (A) 


eps=eps 


toldef=10*max(size(A))*norm(A,1)*eps 
tol=100*eps*norm([A Bl] ) 
[Alc, Ble, Cic,t, k)-ctrbt(A, Bl Cl re) 
tol=100*eps*norm([A B2])) 
[A2c, B2e¢,C2c,t, k)/=cerbt(A, b2 62,05) 
tol=100*eps*norm( [A;C1])) 
[Alo,Blo,Clo,t,k]=obsv£(A,B1,Cl,tol) 
tol=100*eps*norm( [A;C2]) 
[(A20,B20,C20,t,k])=obsvf(A,B2,C2,tol) 


pause 
disp(' 
disp(' 


<) 
") 


P12 


asi ( ‘ ca a a ea a a a a a a a a a a a a a a ee a 


disp(' h2lqg % Running script file H2LQG.M for H2 optimization' ) 

Se a Se eS SS SS = === ") 
aretype='Schur' 

h2lqg 

disp(' ') 

disp(' ") 

disp(' (strike a key to continue ...)') 
clear functions 

pause 

pltopt % Preparing singular values for plotting 

end 

disp(' ) 

disp(' ") 

disp(' (strike a key to continue ...)') 
pause 

disp(' 2) 

Sisp(' ') 

disp(' State space representation of controller (acp, bcp, ccp, dcp) ') 
disp(' anceCiir Tyiul (acl, bel, cel, acl) and poles, controllability, ‘') 
disp(' observability, and condition number of the controller. ') 

< ye) 

disp(' ") 

acp,bcp,ccp,dcp 


polec=eig(acp) 

tol=100*eps*norm([acp bcp] ) 
feaere,bCpc,ccpc,t,k)j=ctrbf(acp,bcp,ccp,tol) 
tol=100*eps*norm( [acp;ccp] ) 
[acpo,bcpo,ccpo,t,k]j=obsvf(acp,bcp,ccp,tol) 
condacp=cond(acp), rcondacp=rcond(acp) 


meme l,ccl,dcl 

polet=eig(acl) 

disp(' ") 

disp(' ") 

disp(' (strike a key to continue ...)') 
pause 

disp(' ") 

disp(' ") 

disp(' Open loop state space representation of controller/plant series') 
disp(' ) 


(algf,blgf,clgf,dlgf]=series(acp,bcp,ccp,dcp,ag,bg,cg,dg) 
polol=eig(algf) 


disp(' ) 

disp(' ") 

disp(' (Serike a Key to Continue ...)  } 
pause 

disp(' ") 

disp(' ') 
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disp(' Closed loop state space representation of controller/plant series,') 


disp(' controllability, observability, and condition number of the closed’) 
disp(' loop acgf matrix. ') 

disp(' ") 

disp(' a) 


[acgf, begf, ccgf,dcgf )=feedbk(algf,blgf,clgf,dlgf,2) 
tol=100*eps*norm({acgf bcgf)) 
{[acgfc,bcgfc,ccgfc,t,k)=ctrbf (acgf,bcgf,ccgf,tol) 
tol=100*eps*norm( {acgf;ccgf] ) 

{acgfo, bcgfo,ccgfo,t,k)=obsvf (acgf,bcgf,ccgf,tol) 
condacgf=cond(acogf ) 


disp(' ") 

disp(' =) 

disp(' (strike a key to continue ...)') 
pause 

disp(' ) 

disp(' Poles of the closed loop system' ) 

disp(' 3) 

polcl=eig(acgf) 

end 
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*x29hinf.m 
diary x29hinf.dat 
format short e 


disp(' ") 

ep(' '‘) 

disp(' This script file is designed to solve the Hinf control ') 

disp(' problem for the X-29. The 14th order FDLTI model, in state ') 
disp(' space form, is that of the X-29 aircraft and actuator dynamics.') 
disp(' Two states are those of the aircraft dynamics, i.e., alpha ') 
disp(' and q. The remaining 12 states are the dynamics of the ') 

disp(' three, fourth order actuators, i.e., the canard, flaperon, and ') 
disp(' strake actuators. The order of the unbalanced states is as ') 
disp(' follows: ') 

disp(' -) 

disp(' alpha, q, dc, df, ds, dcdot, dfdot, dsdot, dcdblidot,') 

disp(' dfdblidot, dsdbldot, dctrpldot, dftrpldot, dstrpldot') 

disp(' wy 

disp(' Given the open loop transfer function G(s)=Cinv(IS-A)B+D, a ') 
disp(' stabilizing controller F(s) will be found such that the Hinf norm’) 
disp(' of Tylul is less than or equal to one. ') 

disp(' ) 

disp(' ") 

pause 

elc 

disp(' X-29 aircraft and actuator state space representation') 

disp(' ") 

disp(' i) 


Sgepeeateola+OO 8 .9960dt00 -.2269d-OlL <-.1213d+00 -.1948d-01 -.9493d-03 .. 
~44270d-04 -.6712d-04 .1451d-05 -.2162d-04 -.3540d-05 0.0 0.0 0.0; 


~-5474d+01 -.3424d+00 Moo satOte— i 3ooa+o! o-.1056datOl .3898d—-Oz2m.. . 
pee sal —.6397d-02 =—.2509d-03 -—.5362d-03 -.2912d-03 0.0 0.0 070; 
mmeeeoc0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 90.0 0.0 0.0 0.0 0.0; 
ameeeo.0 0.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0; 
fee OOO O20) O20" 0.0 0.0 1.0 "0.0 0.0 0.0 0.0 0.0 0.0; 
Zeemeoo oOo O00 0-0 0.2.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0 0.0; 
faemeovono.O) 0-0 0.0 60.0 0.0 0.0 0.0 1.0 0.0 0.0 0.0 0.0; 
fee O-0 0-0 0.0 0.0 0.0 0.0 80.0 0-0) 1.0 0.0%0.0%0.0% 
oe 0-0 0-0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0e+04 0.0 0.0; 
fee o-O 0-0 0.0 0.0 0.0 0.0 0.0 0.0 90.0 0.0 0.0 1.0e+04 0.0; 
peor 0.0 0-0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 1.0e+04; 
moomeo-0 —.14790d+t04 0.0 0.0 -.1143d+03 0.0 0.0 -.2529d+01 0.0... 
Sue —.2097d+03 0.0 0.0; 
mmomeO-0 90.0 —-.1491d+04 0.0 0.0 =-.1149d+03 0.0 0.0 -.2536d+01 ... 
mee 0-0 =-.2701d+03 0.0; 
meeeo-O 0.0 0.0 -—.5302d+05 90.0 0.0 -.1816d+04 0.0 0.0... 


me200tO2 0.0 0.0 -.6053d+03) 


215 


bg={0.0 0.0; 
0.0 Ong 
0.0 0.0; 
O20 “O20, 
0.0 0.0; 
0.0. *0:.. 0; 
020,402.07 
0.0 0.0; 
0-0 020; 
0.0 0.0; 
0.0 0.0; 
-1479d+04 0.0; 
0.0 .1491d+04; 
O00 .£5308a+05) 
eg=(1.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 050 02007) ee 
0.0 1.0 0.0 0.0 0.0 0.0 90.0 0.0 0.0 0.0 0.0 0:0 (0 70RaaeS 
dg=zeros(2) 
pause 
disp(' ") 
disp(' Balanced realization of the X-29 state space representation’ ) 
disp(' ” 
[agbl,bgb1l,cgbl,g,t]=obalreal(ag,bg,cqg) 
ag=agbl; bg=bgbl; cg=cgbl; 
disp(' ") 
disp(' Calculate the poles and transmission zeros of the balanced open’) 
disp(' loop plant’) 
disp(' 4 
poleg=eig(ag), tzerog=tzero(ag,bg,cg,dg) 
disp(' aS 
disp(' “2 
disp(' Determine determine the condition number of ag') 
disp(' “) 
disp(' ") 
condag=cond(ag), rcondag=rcond(ag) 
disp(' 2) 
disp(' -) 
pause 
disp(' ") 


6 


disp(' << Design Specifications >> ') 

disp(' ") 

disp(' 1). Robustness Spec. : -40 dB roll-off, -20 db @ 100 Rad/Sec.') 
disp(' Associated Weighting: ') 

disp(' ") 

disp(' a 1000 ') 

disp(' W300). ————— +r (fixd)') 

disp(' 2 2x2") 

disp(' 8 ") 

disp(' ") 

disp(' ") 

disp(' 2). Performance Spec.: minimizing the sensitivity function') 
disp(' as much as possible.') 

disp(' Associated Weighting: ') 

disp(' i) 

disp(' ) 

disp(' -1 -1 .01(100s + 1) ") 

disp(' W1(s) = Gam * ----------- * T') 

disp(' 2xX2') 

disp(' (.Ols + 1) ") 

disp(' i) 

disp(' where "Gam" in this design is iteratively updated from 1') 
w=logspace(-3,3,100); 

k=1000; mn=[{2 2]; tau=0.0; 

nuw31i = [0.0 k]; dnw3i = [1.0 0 O]; 

SVvw31 bode (nuw3i,dnw31i,w); svw3i = 20*10g10(svw31) ; 

nuwli {1.0 0.01); dnwli =[0.01 1.0); 

Svwli = bode(nuwli,dnwli,w); svwli = 20*10g10(svwli); 

aw2=-0.001*eye(2); bw2=zeros(2); cw2=zeros(2); dw2=-0.001*eye(2); 


disp(' ") 

disp(' ) 

disp(' (strike a key to see the plot of the weightings ...)') 
pause 

semilogx(w,svwli,w, 8vw3i) 

grid 


title('X-29 Design Specifications’ ) 
xlabel('Frequency - Rad/Sec') 
ylabel('1/W1 & 1/W3 - db') 
text(.01,0,'1/W1(s)') 
text(.5,100,‘'1/W3(s)') 

meta x29specs 

pause 


ee) 


Crc 

disp(' 
disp(' we) 
disp(' 
disp(' ) 
disp(' 
disp(' ") 
disp(' 
disp(' i 
disp(' 
disp(' 
disp(' ") 
disp(' 
disp(' 
disp(' —'') 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' i) 
Gita (nas) 
disp(' "ay 
disp(' 
pause 

cle 

disp(' a 
disp(' —_') 
disp(' 
disp(' ") 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 

Gisp (: 
disp(' 
disp(' 
disp(' ") 
disp(' ") 
disp(' 
pause 

exe 

disp(' ") 
disp(' ") 
disp(' =) 


where 


") 


+ + + + + © © F* % & * *& 


<< Problem Formulation >>') 


ye 


2). W3 penalizing plant output "y"') 


F(8) 


Tylul = 


x *« &*& k* *&* &* & & & & & & &K KK kK kK KH KK KK KK KE KK KF 


[Step 1]. 


{Step 2]. 


(Step 3). 
{Step 4]. 


x *& k* & & & kK &e kK &e kK KR Ke Ke KK KK KK KK K Ke KE KE K KK XK 


| Dyas, 


=i 


Gam*W1*(I + GF) 


= 


W3*GF*(I + GF) 


(strike a key to continue 


<< DESIGN PROCEDURE >>') 


Do plant augmentation (run AUGMENT.M or 


AUGX29.M) 


< OF 
inf”) 


"a 


") 


LC cae) 


TS) 


Gam*W1 penalizing error signal "e"') 


and find a stabilizing controller F(s) such that the Hinf 
of TF Tylul is less than or equal to one, 


Gam * Wl * S 


W3 * 


(I - S) 


numerical condition if necessary 


Do Hinf synthesis with "Gam" = 1 


Redo the plant augmentation for a 


higher "Gam" and rerun HINF.M 


(strike a key to continue 


1S 


norm' ) 


- 
ee? ee 


--)') 


Balance the augmented plant for better 


-)9 


Form an augmented plant P(s) with these two weighting functions: ') 


a) 
a 
o 
a) 
>) 
*") 
—_) 
ee 
— 
_ 
ee, 


ary 


disp(' 
disp(' 
disp(' 


Assign the cost coefficients "Gam" with Gam=1 ') 
> 


serving as the baseline design ....') 


gama=[(12.5]; 
ngama=length(gama) ; 
for i=l:ngama 
Gam=gama(1,i) 


k=1000 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
augx29 
disp(' 
disp(' 
disp(' 
pause 
elc 
disp(' 
disp(' 


' ) ' 


(strike a key to continue ...)') 


=) 


tdisp(' 
%disp(' ry 
tdisp(' a) 


Do state space balancing on the augmented plant if needed') 


%{[abal,bbal,cbal,g,t])=obalreal(A, [Bl B2],[(C1;C2)] ) 


%*A=abal, 


disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 


Bl=bbal(:,1:2), B2=bbal(:,3:4), Cl=cbal(1:6,:), C2=cbal(7:8,: ) 

=) 

7) 

The transmission zeros, poles and condition number of the augmented') 
plant follow. In addition, determine if (A,Bl1) & (A,B2) are ') 
stabilizable and if (C1,A) & (C2,A) are detectable.') 

») 

") 


tzeroaug=tzero(A, [Bl B2],[(C1;C2),[D11 D12;D21 D22)), poleaugA=eig(A) 
condaugA=cond(A), rcondaugA=rcond (A) 


eps=eps 


toldef=10*max(size(A))*norm(A,1)*eps 
tol=100*eps*norm([A Bl} ) 
[Aic,Bic,Cic,t,k]=ctrbf(A,B1,C1,tol) 
tol=100*eps*norm([A B2]) 
[A2c,B2c,C2c,t,k)J=ctrbf(A,B2,C2,tol) 
tol=100*eps*norm([A;C1] ) 
[Alo,Blo,Clo,t,kj=obsvf(A,B1,Cl1,tol) 
tol=100*eps*norm( [A;C2]) 
[A20,B20,C20,t,k]=obsvf (A,B2,C2,tol) 


pause 
disp(' 
disp(' 


9) 
=) 


dg 


Gisp(' 0 ter = = = = = - === ") 


disp(' hinf % Running script file HINF.M for Hinf optimization’ ) 
GQisp(* 0 3 9 9 8 = = 8 = eee ") 
aretype='Schur' 

hinf 

disp(' ” 

disp(' ) 

disp(' (strike a key to continue ...)') 
Clear functions 

pause 

disp(' ~) 

disp(' =) 

disp(' State space representation of the full order controller') 
disp(' (acp, bep, ccp, dcp) with its poles and condition number') 
disp(' ") 

acp,bcp,ccp,dcp 


polec=eig(acp) 
condacp=cond(acp) 


disp(' ") 
disp(' Minimal realization of the controller') 
disp(' 2 


toldef=10*max(size(acp))*norm(acp,1)*eps 
tol=100*eps*norm([acp bcp;ccp dcp}) 
(acpm, bcpm, ccpm, dcpmj=minreal(acp,bcp,ccp,dcp) 


disp(' =) 

disp(' Balanced realization & model reduction of the minimal controller') 
disp(' ") 

[acpbl,bcpbl,ccpbl,g,t)=obalreal(acpm,bcpm,ccpm) 

elim=(10) 

(acpr,bcpr,ccpr,dcpr ]=modred(acpbl,bcpbl,ccpbl,dcpm, elim) 

disp(' *) 

disp(‘ Poles, controllability, observability, and condition of the ') 
disp(' balanced, reduced order controller') 

disp(' io 


poleacpr=eig(acpr) 

tol=100*eps*norm([acpr bcpr}) ) 

[acpe,bepce, cepce,t,k])=—ctrb£ (acpr, bepr,cepre, tol) 
tol=100*eps*norm( [acpr;ccpr} ) 

[acpo,bcpo, ccpo,t,k])=obsvf(acpr, bepr cepr, tou) 
condacpr=cond(acpr), rcondacpr=rcond(acpr) 
acp=acpr; bep=bcpr; ccp=ccpr; dcp=dcpr; 

disp(' is 

disp(' CLTF Tylul (acl, bcl, ccl, dcl) and its poles (reduced orderja® 
[acl,bcl,ccl,dcl)=l£ftf (sysp,dimp, acp, bep, cep, acer) 
polet=eig(acl) 


pause 
disp(' zy) 
disp(' oa 
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pltopt % Preparing singular values for plotting 


end 

disp(' ue) 

disp(' ") 

disp(' (strike a key to continue ...)') 
pause 

disp(' = 

disp(' ") 

disp(' Open loop state space representation of controller/plant series’) 
disp(' a) 

[algf,blgf,clgf,dlgf]=series(acp,bcp,ccp,dcp,ag,bg,cg,dg) 

disp(' ") 

disp(' H) 

disp(' (strike a key to continue ...)') 
pause 

disp(' v2) 

disp(' sc) 

disp(' Closed loop state space representation of controller/plant series,') 
disp(' controllability, observability, and condition number of the closed') 
disp(' loop acgf matrix. ') 

disp(' ) 

disp(' ) 


[acgf,bcgf,ccgf,dcgf ]=feedbk(algf,blgf,clgf,dlgf,2) 
tol=100*eps*norm([acgf bcgf]) 
Saegrc,bcofc,ccofc,t,k)]=ctrbf(acgf,bcgf,ccgf,tol) 
tol=100*eps*norm( [acgf;ccgf] ) 
[acgfo,bcgfo,ccgfo,t,k]=obsvf(acgf,bcgf,ccgf,tol) 
condacgf=cond(acgf) 


disp(' ) 

disp(' ") 

disp(' (strike a key to continue ...)') 
pause 

disp(' ") 

disp(' Poles of the closed loop system') 

disp(' ") 

polcl=eig(acgf) 

end 
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% augx29.m 

% Plant Augmentation for the X-29 H2 and Hinf problem as W3 is not a 
% proper transfer function. Includes contingency for adding W2 to 

% ensure dl2 is full column rank. This script file is designed for 
% the X-29 system without theta as a state, ie, 2 inputs & 2 outputs. 
% This script file is a modified version of the aughimat.m M-file taken 
% from the Matlab Robust-Control Tool Box [Ref 7]. 

disp(' ) 

disp(' ") 

disp(' << Plant Augmentation >>') 

Gam=gama(1,i) 

Gam = input(' Input the cost coefficient "Gam" = '); 


cgb = 1/k*[{cg(1,:)*ag*ag;cg(2,:) *ag*ag*ag*tau+cg(2,:)*ag*ag] 

dgb 1/k*[cg(1,:)*ag*bg;cg(2,:)*ag*ag*bg*tau ] 

nwl Gam*{dnw1li;0 0;0 O;dnwli)j 

dQwl = nuwli 

sysw2=[aw2 bw2;cw2 dw2]; xw2=2; 
[A,B1,B2,C1,C2,D11,D12,D21,D22)=augx29pl(ag,bg,cg,cgb,dg,dgb,nwl,dwl, sysw2,xw2 


,mn) 

disp(' ") 

disp(' - - - State-Space (A,B1,B2,C1,C2,D11,D12,D21,D22) is réadyeiaues 
disp(' the Small-Gain problem - - -') 

% 

$ e2---- End of AUGX29.M ----- S 2Z 
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taugx29pl.m 

function 

{a,b1,b2,c1,c2,d11,d12,da21,da22 ]=augx29pl(ag,bg,cga,cgb, dga, dgb, num, den, sys,x,mn) 
% 

% [A,B1,B2,C1,C2,D11,D12,D21,D22 ]=AUGX29PL(AG, BG, CG, CGB, DG, DGB, NUM, DEN, SYS, X, MN) 
% produces the augmented plant for the X-29 h2 and hinf problem for 

an improper W3 weighting and W2 added to make d12 full column rank. 

This script file is a modified version of the augpl.m function taken 
from the Matlab Robust-Control Tool Box [Ref 7]. 


Weighted Plant : 


Xg = Ag Xg + Bg Ug 
iYgaj = iCgaj Xg + {Dga; Ug 
i¥Ygb;  iCgb; iDgb| 
Weighting W1 : Weighting W2 : 
Xwl = Awl Xwl + Bwl Uwl Xw2 = Aw2 Xw2 + Bw2 Uw2 
Ywl = Cwl Xwl + Dwl Uwl Yw2 = Cw2 Xw2 + Dw2 Uw2 


Over all augmented plant : 


cP dP dP dP dP dP dP dP dP dP dh dP dP oP oP dP dP dP dP dP Gg dP cd dP dP dP dP ge oP 


iXg ; | Ag 0 O; iXg i; 10 Bg | ; Ul; 
Tce Tat 1 pXwly + | 1 102; 
'Xwl1! 1-BwlCga Awl Oe XW, 'Bwl -BwlDga| 
cit 
1Xw2 | i Oo 0 Aw2 |! 10 Bw2 
HY deat i-DwlCga Cwl O}; |Xg } 1Dwl -DwlDgaj} ;|U1} 
Ly bb | a 0 Cw2 | {|Xwl} © Dw2 | |U2=Ug} 
wie) = | Cab 0 Of i xXw2ies 10 Dgb ! 
Ye ac! '-Cga 0 oF! a -—Dga |} 

% State Space of Weighting Wl & W2: 

% 

[awl,bwl,cwl,dawl] = tfm2ss(num,den,mn(1),mn(2) ) 

[aw2,bw2,cw2,dw2] = sys2ss(sys,x) 

% 

% State Space of Augmented Plant : 

% 

$$ =----- A matrix : 

% 

[rag,cag] = size(ag); 

[rawl,cawl] = size(awl); 

[raw2,caw2] = size(aw2); 

a = (ag zeros(rag,cawl) zeros(rag,caw2);-bwl*cga awl zeros(rawl,caw2); 


zeros(raw2,Ccag) zeros(raw2,cawl) aw2]); 


12:3 


%$ ----- B matrix : 


[rbg,cbg] = size(bg); 

{[rbwl,cbwl] = size(bwl); 

[ rbw2,cbw2 ] size(bw2); 

bl = [zeros(rbg,cbwl) ;bwl;zeros(rbw2,cbwl1) J; 


b2 = [bg;—-bwl*dga;bw2]; 

% 

a — C matrix : 

% 

[rcgb,ccgb] = size(cgb); 

[rcwl,ccwl] = size(cwl); 

[rcw2,ccw2] = size(cw2); 

cl = [-dwl*cga cwl zeros(rcwl,ccw2);zeros(rcw2,ccgb) zeros(rcw2,ccwl) cw2; 
cgb zeros(rcgb,ccwl) zeros(rcgb,ccw2) J]; 

{rcga,ccga) = size(cga); 

c2 = [-cga zeros(rcga,ccwl) zeros(rcga,ccw2)]; 

% 

% ----- D matrix : 

% 


[rdgb,cdgb] = size(dgb); 

(rdwl,cdw1] = size(dwl); 

(rdw2,cdw2] = size(dw2); 

dll = [dwl;zeros(rdw2,cdwl) ;zeros(rdgb,cdwl) ]; 
dl2 = [-dwl*dga;dw2;dgb]; 


{rdga,cdga] = size(dga); 

d21 = eye(rdga); 

d22 = -dga; 

% 

c= =—--e— End of AUGX29PL.M ---- % 
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$x29analy.m 
diary x29analy.dat 


format 


disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 
disp(' 


short e 


i) 

") 

This script file analyzes the results of Hinf synthesis of the') 
X-29 MIMO model by plotting resultant system return difference’) 
singular values, calculating the closed loop poles/zeros,and ') 
making Bode plots') 

") 

") 

The unaugmented state space plant representation') 

0) 

oy 


ag,bg,cg,dg 


disp(' 
disp(' 
pause 

disp(' 
disp(' 
disp(' 


") 
.) 
State space representation of the controller') 
ia) 
") 


aep, >cp,ccp,dcp 


disp(' 
disp(' 
pause 

disp(' 
disp(' 
disp(' 
disp(' 


r=} 
") 


Open loop state space representation of controller/') 
plant series') 

a 

") 


aor, blgf,clof,dlgf 
[alfg,blfg,clfg,dlfg}=series(ag,bg,cg,dg,acp,bcp,ccp, dcp) 


esp ( * 
disp(' 
disp(' 
disp(' 
disp(' 


7) 
a) 
Compute singular value plot of return difference matrices') 
") 
), 


w=logspace(-3,3,100); 
svg=sigma(ag,bg,cg,dg,3,w); svg=20*1logl0(Ssvg); 
semilogx(w,svg) 

title('X-29 SV PLOT (I + G)') 

xlabel( ‘FREQUENCY - rad/sec') 

ylabel('SV - dB') 


Grad 


meta x29svl 


pause 


svgf=sigma(algf,blgf,clgf,dlgf,3,w); svgf=20*logl0(svgf); 
semilogx(w,svgf ) 

title('X-29 SV PLOT (I + GF)') 

xlabel('FREQUENCY - rad/sec') 

ylabel('SV - dB') 


gira 


meta x29sv2 
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pause 
svfg=sigma(alfg,blfg,clfg,dlfg,3,w); svfg=20*logl0(svfg); 
semilogx(w, svfg) 

title('X-29 SV PLOT (I + FG)') 

xlabel('FREQUENCY - rad/sec') 

ylabel('SV - dB") 

Gis 

meta x29sv3 

pause 

svgf=sigma(algf,blgf,clgf,dlgf,4,w); svgf=20*logl0(svgf) ; 
semilogx(w, svgf) 

title('X-29 SV PLOT (I + inv(GF))') 

xlabel('FREQUENCY - rad/sec') 

ylabel('SV - dB') 

grid 

meta x29sv4 

pause 

svfg=sigma(alfg,blfg,clfg,dlfg,4,w); svfg=20*logl0O(svfg); 
semilogx(w,svfg) 

title('X-29 SV PLOT (I + inv(FG))') 

xlabel('FREQUENCY - rad/sec') 

ylabel('SV —- dB') 

grid 

meta x29sv5 

pause 

svloop=sigma(algf,blgf,clgf,dlgf,1,w); svloop=20*log10(svloop) ; 
semilogx(w,svloop) 

title('X-29 SV PLOT (GF)') 

xlabel('FREQUENCY - rad/sec') 

ylabel('SV - dB') 


grid 

meta x29sv6 

pause 

disp(' ") 

disp(' ) 

disp(' Closed loop state space representation of controller/plant series’) 
disp(' ") 

disp(' ") 

acgft, Degft,ccegt jacgt 
tzerocgf=tzero(acgf,bcgf,ccgf,dcgf) 
disp(' ") 

disp(' ") 

pause 
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disp(' Poles and zeros of the closed loop controller/plant system ') 


disp(' Waegt, begf, ccgf, dgef) *) 
disp(' =) 

disp(' =) 

disp(' output alpha/q vs input 1 ') 
disp(' ") 
{z,p,k]=ss2zp(acgf,bcgf,ccgf,dcgf,1) 
pause 

disp(' “4 

disp(' output alpha/gq vs input 2 ') 
disp(' +) 
{(z,p,k])=ss2zp(acgf,bcgf,ccgf,dcgf,2) 
pause 

disp(' ") 

disp(' ") 

disp(' Open loop Bode plots of outputs vs inputs') 
disp(' :s) 

disp(' ") 


{magl, phasel]=bode(algf,blgf,clgf,dlgf,1,w); mag1l=20*1l0g10(mag1l) ; 
{mag2,phase2 ]=bode(algf,blgf,clgf,dlgf,2,w); mag2=20*1l0g10(magZ2) ; 
semilogx(w,magl(:,1)) 

title('X-29 OPEN LOOP BODE PLOT INPUT 1 / alpha (GF)') 
xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 

grid 

meta x29bodel 

pause 

semilogx(w,magl(:,2)) 

title('X-29 OPEN LOOP BODE PLOT INPUT 1 /q _ (GF)') 
xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 

quar 

meta x29bode2 

pause 

semilogx(w,mag2(:,1)) 

title('X-29 OPEN LOOP BODE PLOT INPUT 2 / alpha (GF)') 
xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 

grid 

meta x29bode3 

pause 

semilogx(w,mag2(:,2) ) 

title('X-29 OPEN LOOP BODE PLOT INPUT 2 / q_ (GF)") 
xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 


grid 

meta x29bode4 
disp(' ap 
disp(' *) 
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disp(' Closed loop Bode plots of outputs vs inputs’) 
disp(' ") 

disp(' ") 

[mag1, phasel )=bode(acgf,bcgf,ccgf,dcgf,1,w); magl=20*10g10(magl1); 
[mag2, phase2 ]=bode(acgf,bcgf,ccgf,dcgf,2,w); mag2=20*1l0g10(mag2) ; 
semilogx(w,magl(:,1)) 

title('X-29 CLOSED LOOP BODE PLOT INPUT 1 / alpha ‘') 
Xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB") 

grid 

meta x29bode5 

pause 

semilogx(w,magl(:,2)) 

title('X-29 CLOSED LOOP BODE PLOT INPUT 1/q ') 
xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB’) 

grid 

meta x29bode6 

pause 

semilogx(w,mag2(:,1)) 

title('X-29 CLOSED LOOP BODE PLOT INPUT 2 / alpha ‘') 
xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 

grid 

meta x29bode7 

pause 

semilogx(w,mag2(:,2)) 

title('X-29 CLOSED LOOP BODE PLOT INPUT 2 /q ') 
xlabel('FREQUENCY - rad/sec') 

ylabel('GAIN - dB') 

grid 

meta x29bode8 

pause 

end 
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%x29resp.m 
diary x29resp.dat 
format short e 


disp(' ") 

disp(' ) 

disp(' This script file is designed to calculate and plot the time') 
disp(' response of the augmented X-29 (controller/plant closed loop') 
disp(' series) to a 1 degree / 1 second step input from each of') 
disp(' the two inputs separately.') 

disp(' =} 

disp(' ) 

disp(' Setting up the time vectors') 

disp(' nt) 


time=(0:0.01:4] 

tstas=(10000;0 1000;0 0100;000 10;0000 1) 
*ccgfwc=[zeros(5,18) stas zeros(5,6) ] 

tdcgfwc=zeros(5,2) 


stas=(10000000;01000000;00100000;0001000 Q;... 
MeemoeO 1 000;00000100;000000107;000000 0 1) 
ccgfwc=[zeros(8,15) stas zeros(8,6) ] 

dcgfwc=zeros (8,2) 

disp(' z) 

disp(' Plotting alpha, q, dc, df, ds, dcdot, dfdot, and dsdot') 

disp(' responses to a 0.01745 rad / 1 second step from input one') 
disp(' ‘)) 


u={0.01745*ones(1,101) zeros(1,300);zeros(1,401)]' 
{yJ=lsim(acgf,bcgf,ccgfwc,dcgfwc,u,time) 
${yJ=lsim(acgf,bcgf,ccgf,dcgf,u,time) 

plot(time,y(:,1) ) 

title('X-29 alpha RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 1)') 
xlabel('TIME - sec') 

ylabel('radians' ) 

grid 

meta x29rspll 

pause 

plot(time,y(:,2) ) 

title('X-29 q RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 1)') 
xlabel('TIME - sec') 

ylabel('radians/second' ) 

grid 

meta x29rspl2 

pause 

plot(time,y(:,3) ) 

title('X-29 DC FOR 0.01745 rad / 1 sec STEP (INPUT 1)') 
xlabel('TIME - sec') 

ylabel('radians' ) 

grid 

meta x29rsp13 

pause 
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plot (time,y(:,4)) 

title('X-29 DF FOR 0.01745 rad / 1 sec STEP (INPUT 1)"') 
xlabel('TIME - sec') 

ylabel(‘radians' ) 

grid 

meta x29rspl14 

pause 

plot(time,y(:,5)) 

title('xX-29 DS FOR 0.01745 rad / 1 sec STEP (INPUT 1)') 
xlabel('TIME - sec') 

ylabel('‘radians') 

grid 

meta x29rsp15 

pause 

plot(time,y(:,6) ) 

title('X-29 DCDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1)') 
xlabel('TIME - sec") 

ylabel('radians/second' ) 

grid 

meta x29rspl16 

pause 

plot (time,y(:,7)) 

title('X-29 DFDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1)') 
xlabel('TIME - sec') 

ylabel('radians/second') 

grid 

meta x29rsp17 

pause 

plot (time;y (=>) 

title('X-29 DSDOT FOR 0.01745 rad / 1 sec STEP (INPUT 1)') 
xlabel('TIME - sec') 

ylabel('radians/second' ) 


grid 

meta x29rspl18 

pause 

disp(' ") 

disp(' Plotting alpha, q, dc, df, ds, dcdot, dfdot, and dsdot') 
disp(' responses to a 0.01745 rad / 1 second step from input two') 
disp(' ") 


u=(zeros(1,401);0.01745*ones(1,101) zeros(1,300) ]' 
[yJ=lsim(acgf,bcgf,ccgfwc,dcgfwc,u,time) 
*{yJ=lsim(acgf,bcgf,ccgf,dcgf,u,time) 

plot(time,y(:,1)) 

title('X-29 alpha RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 2)') 
xlabel('TIME - sec') 

ylabel('radians') 

grid 

meta x29rsp21 

pause 
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plot (time,y(:,2)) 

title('X-29 q RESPONSE TO 0.01745 rad / 1 sec STEP (INPUT 2)') 
xlabel( ‘TIME - sec‘) 

ylabel(‘'radians/second' ) 

chante 

meta x29rsp22 

pause 

plot (time,y(:,3)) 

title('x-29 DC FOR 0.01745 rad / 1 sec STEP (INPUT 2)') 
xlabel('TIME - sec') 

ylabel('radians') 

grid 

meta x29rsp23 

pause 

plot(time,y(:,4)) 

title('X-29 DF FOR 0.01745 rad / 1 sec STEP (INPUT 2)') 
xlabel('TIME - sec') 

ylabel('radians' ) 

grid 

meta x29rsp24 

pause 

plot(time,y(:,5)) 

title('X-29 DS FOR 0.01745 rad / 1 sec STEP (INPUT 2)') 
xlabel('TIME - sec') 

ylabel(‘radians' ) 

grid 

meta x29rsp25 

pause 

plot (time,y(:,6)) 

title('X-29 DCDOT FOR 0.01745 rad / 1 sec STEP (INPUT 2)'") 
xlabel('TIME - sec") 

ylabel('radians/second' ) 

grid 

meta x29rsp26 

pause 

plot(time,y(:,7) ) 

title('X-29 DFDOT FOR 0.01745 rad / 1 sec STEP (INPUT 2)') 
xlabel('TIME - sec') 

ylabel('radians/second' ) 

grid 

meta x29rsp27 

pause 

plot (time,y(:,8)) 

title('X-29 DSDOT FOR 0.01745 rad / 1 sec STEP (INPUT 2)') 
xlabel('TIME - sec") 

ylabel('radians/second' ) 

grid 

meta x29rsp28 

end 
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%*x29schmr.m 
Giary x29schmr.dat 
format short e 


disp(' ") 

disp(' a) 

disp(' This script file employs Schur Model Reduction ') 

disp(' to reduce the order of the minimal, balanced hinf controller. ') 
disp(' The reduced order controller is then combined in series with the ') 
disp(' balanced X-29 plant. The reduced order model and minimal, balanced') 
disp(' model SV Bode plots are compared and the error system plotted.') 
disp(' Additionally, the reduced order controller/balanced plant series’) 
disp(' is compared with the minimal, balanced controller/plant series.') 
disp(' ") 

disp(' ) 

disp(' State space representation of the minimal, balanced hinf controller' ) 
disp(' ") 

acp,bcp,ccp,dcp 


poleacp=eig(acp) 
condacp=cond (acp) 


disp(' ') 
disp(' ") 
adisp(' Schur Model Reduction of the hinf controller') 
disp(' ") 


{[acph, bcph, ccph,dcph,totbnd, hsv)=schmr(acp,bcp,ccp,dcp,1,10) 
poleacph=eig(acph) 
condacph=cond(acph) 


disp(' ") 
disp(' Computing the SV Bode plot of the minimal & reduced order controllers’ ) 
disp(' *) 


w=logspace(-3,3,100); 

svCp=sigma(acp,bcp,ccp,dcp,1l,w); svcp=20*1logl0(svcp); 
svcph=sigma(acph,bcph,ccph,dcph,1l,w); svcph=20*1log10(svcph); 
{acpher,bcpher,ccpher,dcpher ])=addss(acp,bcp,ccp,dcp,acph,bcph,-ccph,-dcph) ; 
sver=sigma(acpher,bcpher,ccpher,dcpher,1,w); sver=20*logl0(sver); 
totbnd=20*1l0g10(totbnd) *ones(w); 

semilogx(w,svcp) 

title('15th ORDER X-29 CONTROLLER’ ) 

xlabel('FREQUENCY - rad/sec') 

ylabel('SV - dB') 

grad 

meta x29schrl 

pause 

semilogx(w, s8vcp,w,S8vcph) 

title('10th & 15th ORDER X-29 CONTROLLERS' ) 

xlabel('FREQUENCY - rad/sec') 

ylabel('SV - dB') 

grid 

meta x29schr2 

pause 
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semilogx(w,totbnd,w,sver) 

title('10th ORDER X-29 CONTROLLER ERROR & ERROR BOUND') 
xlabel( ‘FREQUENCY - rad/sec') 

ylabel('SV - dB') 

text(.02,6,'error bound' ) 


grid 

meta x29schr3 

pause 

disp(' ") 

disp(' State space representation of the minimal, balanced controller/plant') 
disp(' series and ROM controller/balanced plant series (open loop)') 

disp(’ ") 

algf,blgf,clgf,dlgf 


polealgf=eig(algf) 

condalgf=cond(algf) 
[algfrc,blgfrc,clgfrc,dlgfrc]=series(acph,bcph,ccph,dcph,ag,bg,cg,dg) 
polealgfrc=eig(algfrc) 

condalgfrc=cond(algfrc) 


disp(' ") 

disp(' Computing the SV Bode plot of the minimal, balanced controller/plant') 
disp(' series and the ROM controller/balanced plant series.') 

disp(' ") 


svlgf=sigma(algf,blgf,clgf,dlgf,1,w); svlgf=20*loglO(svlgf); 
svrc=sigma(algfrc,blgfrc,clgfrc,dlgfrc,1,w); svrc=20*loglO(svrc); 
{[algfer,blgfer,clgfer,dlgfer]j=addss(algf,blgf,clgf,dlgf,algfrc,blgfrc,... 
eeamnnc,—-Algfrc) ; 

svlgfer=sigma(algfer,blgfer,clgfer,dlgfer,1,w); svlgfer=20*logl0(svlgfer) ; 
semilogx(w,svligf,w,svrc) 

title('1l0th & 15th ORDER X-29 CONTROLLERS IN SERIES WITH X-29 PLANT') 
xlabel('FREQUENCY - rad/sec') 

ylabel('SV - dB') 

grid 

meta x29schr4 

pause 

semilogx(w,svlgfer) 

title('1l0th ORDER X-29 CONTROLLER/PLANT SERIES ERROR’) 

xlabel( ‘FREQUENCY - rad/sec') 

ylabel('SV - dB') 

grid 

meta x29schr5 

end 
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APPENDIX B 


X-29 MODEL STATE SPACE REALIZATIONS 


Open loop state space representation of the uncompensated 


X-29 model. 

ag = 
Columns 1 through 6 
-4.1810e-01 Se Jocoe-01 


5.4740e+00 -3.4240e-01 
0 


O 


© .© 0O7@ © 2 C2 OGeq © © 
Oo 0 0 0 @ 0 OOO © © 


Columns 7 through 12 


4.4270e-05 -6.7120e-05 
=~l.1640e-02 > -6-3970e-03 


0 0 

1 .0000e+00 0 
0 1.0000e+00 

0 0 

0 ) 

0 0 

0 0 

0 0 

0 0 

0 0 
-1.1490e+02 O 


O| 6 =1.8160eE703 


-2.2690e-02 <-1.2130e-01 -1.9480e-02 
2.5850e+00 -1.3860e+00 -1.0580e+00 
O 


© 


) 
0 
0 
0 
0 
0 
0 
0 
0 
0 


-1.4790e+0 
-1.4910e+0 


oO wo OF OC O7e70 6 


O 
O 
O 
O 
0 
0 
0 
0 
0 
3 O 
0 ~ 


=5. 90Z20eC40 


1.4510e-06 -2.1620e-05 -3.5400e-06 
-2.5090e-04 -5.3620e-04 -2.9120e-04 


0 O 0 

0 0 0 

0 0 0 
1.0000e+00 0 O 
O 1.0000e+00 0 

0 0 1.0000e+00 

0 O 0 

0 0 0 

0 0 0 
-2.5290e+00 0 0 
O -2.5360e+00 0 

0 O -1.7900e+01 


134 


~9.4930e=G28 
3.8980e-03 
1.0000e+00 
0 


-1.1430e+0 


OVO. Oss ©. O17 


Ov® OO 0.0 0 Oo 


1.0000e+04 
O 

0 
-2.6970e+02 
O 

O 


bg 


cg 


Columns 13 through 14 


0 0 

0 0 

0 0 

0 0 

0 0 

0 0 

@) 0 

0 0 

0 0 
1.0000e+04 0 
0 1.0000e+04 

0 0 
-2.7010e+02 0 


1.4790e+03 


0 -6.0530e+02 


oo oo o Qo 0 000 O 
oo oo @ © © fo OC fo 2 


© 


0 1.4910e+03 
0 5.3080e+04 


Columns 1 through 12 


i 
0 


Columns 


0 
0 


0 0 0 0 0 
1 0 0 0 8) 


13 through 14 


0 
0 


3.5 


Aa eo, 


Closed loop state space representation of the H, compensated 


acgf = 


Columns 1 through 6 


-1.0380e-02 
-7.4356e-04 
-9.6973e-02 
6.6210e-01 
8.7596e-01 
-1.0283e-01 
1.6027e-01 
-1.1326e-01 
2.-0666e-02 
-2.9553e-02 
6.6863e-03 
2.6231e-03 
6.2838e-03 
-2.3165e-03 
-2.9225e-04 
1.9047e+00 
1.6249e+00 
-5.9083e+00 
-4.0509e+00 
-4.2915e+00 
5.5297e+00 
-2.4640e+00 
-9.0592e-01 
1.8418e+00 
-2.6872e+00 
4.9869e-01 
2.7973e-01 
-1.8133e-01 
3.7460e-02 


2.5396e-04 
-1.0502e-02 
1.5535e-01 
1.2903e+00 
1.1040e+00 
2.-6001le-01 
3.0656e-01 
7.9258e-02 
8.941le-02 
2.4729e-02 
3.3912e-02 
3.0885e-03 
4.2536e-03 
1.5649e-03 
-8.3162e-04 
8.3656e-01 
-6.9125e+00 
7.4810e+00 
3.3677e+00 
-3.9162e-01 
3.0377e+00 
7.1488e-01 
-1.0399e+00 
4.5560e-01 
-7.4178e-01 
1.6324e-01 
1.2875e-01 
9.1365e-02 
1.4831e-02 


-2.0139e-02 
-3.6887e-01 
-2.5920e+00 
1.1697e+02 
8.0559e+01 
-1.9307e+00 
9.7828e+00 
-5.0472e+00 
1.4269e+00 
-1.2592e+00 
4.9473e-01 
1.3971e-01 
3.1809e-01 
-1.0068e-01 
-1.8149e-02 
1.3368e+01 
-8.4058e+01 
8.4663e+01 
3.5998e+01 
-1.1427e+01 
4.6432e+01 
5.2026e+00 
-1.4394e+01 
8.5033e+00 
-1.3371le+01 
2.8015e+00 
2.0369e+00 
8.6801le-01 
2.4260e-01 


1.e., optimal-performance model. 


3.3550e-01 
-2.2137e+00 
-1.1063e+02 
-1.5192e+02 
-2.7211e+02 
8.7694e+01 
-6.1279e+01 
5.0715e+01 
-9.0525e+00 
1.311 76701 
-2.8382e+00 
-1.1509e+00 
-2.7681e+00 
1.0309e+00 
1.2658e-01 
2.6068e+01 
-5.2856e+02 
6.4688e+02 
3.1630e+02 
4.9115e+01 
1.1967e+02 
9.6072e+01 
-5.8765e+01 
-3.1257e-01 
-5.1094e+00 
2.7956e+00 
4.2538e+00 
9.8691e+00 
3.9552e-01 
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6.9246e-01 
-1.9108e+00 
-6.1722e+01 
-1.7331e+02 
-3.5450e+02 
1.8904e+02 
-1.1054e+02 
8.8977e+01 
-1.7498e+01 
2.22176+01 
-5.6200e+00 
-2.0618e+00 
-4.8951e+00 
1.7591e+00 
2.3676e-01 
-1.0156e+01 
-5.2865e+02 
7.1853e+02 
3.7254e+02 
1.2470e+02 
1.2039e+01 
1.3567e+02 
-3.8942e+01 
-3.3913e+01 
4.4225e+01 
-6.4631e+00 
-1.0900e+00 
1.2627e+01 
-3.1035e-01 


-7.5674e-02 
-3.3622e-01 
-4.9836e+00 
-6.7749e+01 
-1.4092e+02 
-8.6324e+00 
1.4974e+01 
-1.7586e+01 
-1.1775e-01 
-4.2756e+00 
-1.2383e-01 
2.3346e-01 
6.4728e-01 
-3.2588e-01 
-1.2454e-02 
1.7859e+01 
-7.8929e+01 
6.9016e+01 
2.5569e+01 
-2.1801l1e+01 
5.9368e+01 
-9.1372e-01 
-1.6421e+01 
1.2907e+01 
-1.9782e+01 
3.9870e+00 
2.6956e+00 
4.1132e-01 
3.3121e-01 


Columns 7 through 12 


7.3362e-02 
-4.7184e-01 
-7.9436e+00 
-7.2469e+01 
-1.0508e+02 
-1.4401e+01 
-8.7251e+01 
8.7848e+01 
-3.6523e+01 
1.8198e+01 
-1.0989e+01 
-2.6019e+00 
-5.6900e+00 
1.5521e+00 
3.7429e-01 
1.1270e+01 
-1.5310e+02 
1.8002e+02 
8.5846e+01 
6.4678e+00 
4.5711le+01 
2.3763e+01 
-1.9057e+01 
3.3588e+00 
-6.5445e+00 
1.7603e+00 
1.7807e+00 
2.-5756e+00 
1.8703e-01 


1.0724e-01 
1.395le-01 
4.4151e+00 
-8.9555e+00 
-4.8144e+01 
1.6726e+01 
-3.1218e+01 
-8.6275e+01 
2.3802e+01 
-4.9880e+01 
2.5817e+00 
3.0928e+00 
7.9380e+00 
-3.5151e+00 
-2.5131le-01 
-1.5754e+01 
3.9986e+01 
-2.1721le+01 
-2.5523e+00 
2.5033e+01 
-5.0002e+01 
7.7899e+00 
1.1990e+01 
-1.2758e+01 
1.9154e+01 
-3.7337e+00 
-2.3549e+00 
3.0174e-01 
-2.9846e-01 


5.6964e-03 
9.9626e-02 
1.8861e+00 
1.2570e+01 
1.3934e+01 
4.4932e+00 
3.5283e+01 
-3.9404e+01 
-2.8890e+01 
1.1589e+01 
-2.3453e+01 
-3.2759e+00 
-6.4217e+00 
8.8043e-01 
6.0387e-01 
-4.6012e+00 
3.2641le+01 
-3.4040e+01 
-1.4893e+01 
3.2072e+00 
-1.6278e+01 
-2.6644e+00 
5.2666e+00 
-2.7551e+00 
4.3892e+00 
-9.3719e-01 
-7.0399e-01 
-3.819le-Ol 
-8.2715e-02 


-3.1043e-02 
2.1029e-02 
-9.1523e-02 
1.0305e+01 
2.2625e+01 
-1.8884e+00 
1.7517e+01 
3.9886e+01 
-7.8925e+00 
-5.9065e+01 
5.7327e+00 
7.3988e+00 
2.0670e+01 
-9.0745e+00 
-6.0252e-01 
2.0579e+00 
8.8556e+00 
-1.5764e+01 
-9.1686e+00 
-6.0269e+00 
5.4075e+00 
-4.3353e+00 
-3.81lle-Ol1 
2.3189e+00 
-3.3116e+00 
5§.9074e-01 
2.9675e-01 
-3.5512e-01 
4.1984e-02 
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-2.5046e-03 
-3.6276e-02 
-6.9471e-01 
-4.3749e+00 
-4.6125e+00 
-1.6597e+00 
-1.0706e+01 
1.2111e+01 
2.3488e+01 
8.0195e-01 
-3.5512e+01 
-1.2501e+01 
-2.1244e+01 
2.0510e+00 
2.0077e+00 
1.7023e+00 
-1.1758e+01 
1.2173e+01 
5.2950e+00 
-1.2489e+00 
5.9970e+00 
9.1071le-01 
-1.9217e+00 
1.0341e+00 
-1.6422e+00 
3.4906e-01 
2.6021e-01 
1.3415e-01 
3.0670e-02 


2.2455e-03 
4.2633e-03 
1.1167e-01 
4.5774e-02 
-6.6097e-01 
3.6889e-01 
5.6280e-01 
-3.9614e+00 
-2.9093e+00 
6.3343e+00 
1.2254e+01 
-9.3120e+00 
-6.7622e+01 
1.3813e+01 
2.4844e+00 
-3.8322e-01 
1.2624e+00 
-9.1123e-01 
-2.5765e-01 
5.5221le-01 
-1.2395e+00 
1.2121le-01 
3.1606e-01 
-2.9692e-01 
4.4927e-01 
-8.8705e-02 
-5.7508e-02 
8.4243e-04 
-7.1986e-03 


Columns 13 through 18 


-5.350le-03 
9.4699e-03 
8.8814e-02 
2.5904e+00 
4.9202e+00 

-1.0100e-01 
5 .0603e+00 
4.5223e+00 

-3.8990e+00 

-2.2728e+01 
1.4122e+01 
4.9965e+01 

-9.5020e+01 
6.1890e+01 
9.9712e+00 
1.1827e-01 
3.4434e+00 

-4.7831e+00 

-2.5074e+00 

-9.2095e-01 
7.6436e-02 

-9.4065e-01 
2.2512e-01 
2.6923e-01 

-3.5903e-01 
5.5418e-02 
1.4788e-02 

-8.6209e-02 
3.0370e-03 


2.0318e-03 
3.5463e-03 
9.5436e-02 
-1.6803e-03 
-6.5175e-01 
3.2152e-01 
4.0275e-01 
-3.4922e+00 
-2.5071e+00 
6.1784e+00 
8.7487e+00 
-1.7175e+01 
2.1892e+01 
-1.2681le+02 
-1.9078e+01 
-3.3412e-01 
1.0399e+00 
-7.1415e-01 
-1.836le-01 
4.9337e-01 
-1.0758e+00 
1.2001le-01 
2.7045e-01 
-2.6170e-01 
3.9521le-01 
-7.7786e-02 
-5.0093e-02 
2.0978e-03 
-6.2893e-03 


1.7690e-04 
-9.9307e-04 
-1.5234e-02 
-1.791le-01 
-2.7845e-01 
-2.3626e-02 
-3.8818e-01 
1.8122e-02 
5.1219e-01 
9.5196e-01 
-1.6967e+00 
5.0142e-01 
1.0389e+01 
6.6685e+00 
-1.4357e+02 
2.3620e-02 
-3.3733e-01 
3.9905e-01 
1.9103e-01 
1.6779e-02 
9.7119e-02 
5.3684e-02 
-4.1326e-02 
6.2769e-03 
-1.2770e-02 
3.5689e-03 
3.7449e-03 
5.7673e-03 
3.8849e-04 


-9.5355e+00 
7.6335e-02 
-4.3986e+01 
3.1638e+02 
4.1820e+02 
-5.2255e+01 
6.9366e+01 
-5.3946e+01 
7.8527e+00 
-1.4194e+01 
2.4494e+00 
1.1575e+00 
2.8333e+00 
-1.1049e+00 
-1.1957e-01 
-3.6424e-02 
-1.5021e+00 
-1.0712e+00 
-1.9232e-01 
2.9619e-01 
-7.1159e-01 
2.2932e-02 
1.8364e-01 
-1.5060e-01 
2.4407e-01 
-5.2288e-02 
-3.1193e-02 
-2.4991le-03 
-4.6699e-03 
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3.4391e+00 
7.7386e+00 
4.4315e+01 
-1.0033e+02 
-1.9565e+02 
-4.8502e+01 
-1.1341le+02 
-2.1608e+01 
-3.4635e+01 
-6.7545e+00 
-1.2916e+01 
-1.3352e+00 
-2.0479e+00 
-3. 92:1 8e-—02 
3.2707e-01 
-1.9965e+00 
-1.7400e+00 
2.2889e+00 
8.9650e-01 
-3.5656e-01 
1.2877e+00 
1.1616e-01 
-3.8875e-01 
2.4093e-01 
-3..8251e=08: 
7.8130e-02 
5.7132e-02 
2.0851le-02 
6.7562e-03 


4.5069e+00 

7.6895e+00 

4.9092e+01 
-1.3583e+02 
-2.4224e+02 
-4.2299e+01 
-1.2071e+02 
-1.5353e+01 
-3.5348e+01 
-5.1031e+00 
-1.3127e+01 
-1.4601e+00 
-2.3598e+00 
-2.6334e-01 

3.3898e-01 

1.5995e-01 
-2.6733e+00 
-1.4713e+0l 
-1.5627e+01 
-6.6464e-01 
-8.3136e+00 
-3.8392e+00 

3.4413e+00 
-8.2107e-01 

1.4705e+00 
-3.6707e-01 
-3.4273e-01 
-4.2393e-01 
-3.690le-02 


Columns 19 through 24 


-2.8169e+00 
-3.5970e+00 
-2.6254e+01 
8.7045e+01 
1.4443e+02 
1.5952e+01 
6.1686e+01 
3.2024e+00 
1.7142e+01 
1.3408e+00 
6.3315e+00 
7.6968e-01 
1.3152e+00 
4.1646e-02 
-1.6766e-01 
-7.3858e-02 
1.2795e+00 
1.5508e+01 
-4.5349e+01 
-1.2548e+00 
-7.4415e+01 
-2.2661le+01 
1.7982e+01 
-2.8087e+00 
6.0983e+00 
-1.7819e+00 
-1.8477e+00 
-2.7233e+00 
-1.9364e-01 


-4.4292e+00 
-8.8914e-02 
-2.0887e+01 
1.4674e+02 
1.9497e+02 
-~2.3194e+01 
3.3636e+01 
-2.4400e+01 
4.1569e+00 
-6.4030e+00 
1.3305e+00 
5.5234e-01 
1.3325e+00 
-5.0056e-01 
-6.0088e-02 
5.4248e-02 
4.3523e-01 
9.3482e-01 
-9.8406e+00 
-1.5613e+01 
5.5042e+01 
-1.3867e+01 
-1.1310e+01 
1.5640e+01 
-2.2100e+01 
4.1876e+00 
2.5316e+00 
-8.2321e-01 
3.2615e-01 


2.5190e+00 

3.6727e+00 

2.5148e+01 
-7.7029e+01 
-1.3179e+02 
-1.8221le+01 
-6.0354e+01 
-5.2758e+00 
-1.7197e+01 
-1.8965e+00 
-6.3686e+00 
-7.4218e-01 
-1.2364e+00 
-~8.3625e-02 

1.666le-01 

8.0145e-02 
-1.2763e+00 
-1.4353e+01 

1.0983e+02 
-7.2403e+00 
-1.8958e+02 
-6.3472e+01 

1.6003e+02 
-9.0324e+01 

1.3275e+02 
-2.5282e+01 
-1.7125e+01 
-3.6217e+00 
-2.0869e+00 


-2.6330e-01 
-1.1671e+00 
-5.4978e+00 
6.6623e+00 
1.8295e+01 
8.6966e+00 
1.5222e+01 
4.6927e+00 
5.0048e+00 
1.3957e+00 
1.8791e+00 
1.7013e-01 
2.3270e-01 
8.8378e-02 
-4.6050e-02 
-3.2383e-02 
3.5565e-01 
4.5759e+00 
-2.5846e+01 
1.2733e+00 
1.4305e+02 
-4.5635e+01 
-1.0863e+02 
4.0505e+01 
-4.7354e+01 
4.8131e+00 
5.4035e-01 
-8.5471e+00 
1.8733e-01 


eS 


1.5369e-01 
1.1007e+00 
4.7458e+00 
-3.1450e+00 
-1.3100e+01 
-8.7138e+00 
-1.3659e+01 
-4.9571e+00 
-4.6390e+00 
-~1.4560e+00 
-1.7467e+00 
-1.4888e-01 
-1.9124e-01 
-9.4238e-02 
4.2216e-02 
3.1757e-02 
-3.2666e-01 
-4.2979e+00 
2.4917e+01 
1.7900e+00 
-1.6158e+02 
1.7203e+02 
-5.1421e+01 
5§.3125e+01 
-7.4804e+01 
1.9246e+01 
1.4710e+01 
5.9972e+00 
1.7541e+00 


-1.9085e+00 
-4.9087e-01 
-1.0658e+01 
6.2425e+01 
8.6622e+01 
-6.0691e+00 
1.9644e+01 
-8.1206e+00 
3.6444e+00 
-2.0670e+00 
1.2745e+00 
2.9148e-01 
6.3395e-01 
-1.6967e-0O1 
-4.2438e-02 
9.4644e-03 
3.1595e-01 
2.1484e+00 
-1.5217e+01 
-1.4103e+01 
7.4043e+01 
-5.2529e+01 
1.9684e+01 
-6.7007e+01 
1.7661e+02 
-3.8141le+01 
-2.4872e+01 
1.1996e+00 
-3.2227e+00 


Columns 25 through 29 


7.6745e-01 
1.2606e+00 
8.1808e+00 
-2.3217e+01 
-4.0969e+01 
-6.7799e+00 
-2.0000e+01 
-2.3565e+00 
-5.8195e+00 
-7.9439e-01 
-2.1598e+00 
-2.4286e-01 
-3.9539e-01 
-3.9883e-02 
5.5939e-02 
2.8776e-02 
-4.2890e-01 
-4.9759e+00 
3.1196e+01 
5.1872e+00 
-1.3929e+02 
1.0694e+02 
-1.0021e+02 
5.1785e+01 
-3.0613e+02 
1.0064e+02 
7.8071e+01 
-3.2832e-01 
9.8153e+00 


3.9425e-03 
=-2.5128e-01 
-9.0223e-01 
-5.7644e-01 
1.2770e+00 
2.-2004e+00 
2.8308e+00 
1.3508e+00 
1.0256e+00 
3.9001le-01 
g-cozoe-OrT 
2.92Z213e-02 
3.2021e-02 
2-6004e-02 
=" 1 300eC-03 
-7.7507e-03 
7.0974e-02 
9.7110e-01 
-5.9282e+00 
1.8317e-01 
2.9207e+01 
-1.9838e+01 
2.2620e+01 
-1.6217e+01 
9.6944e+01 
-5.7549e+01 
-8.9044e+01 
-5.9329e+00 
-9.4011e+00 


-4.4244e-01 
-2.1089e-03 
-2.0616e+00 
1.4670e+01 
1.9437e+01 
-2.3756e+00 
3.2828e+00 
-2.4731e+00 
3.8750e-01 
-6.4995e-01 
1.2241le-01 
5.4373e-02 
1.3221le-01 
-5.0690e-02 
-5.7545e-03 
5.6167e-03 
4.1538e-02 
6.5646e-02 
-9.0624e-01 
-3.1417e+00 
3.1780e+00 
-1.1774e-01 
-6.1728e-01 
-2.1071le+01 
6.6540e+00 
3.7812e+01 
-7.2702e+01 
-4.0260e+00 
-1.7763e+01 
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-2.4640e-01 
5 .5342e-02 
-9.4110e-01 
8.2701e+00 
1.0498e+01 
-1.8131e+00 
1.1850e+00 
-1.6762e+00 
-1.5636e-02 
-4.4838e-01 
-~1.9394e-02 
2.3603e-02 
6.6164e-02 
-3.3977e-02 
-1.1385e-03 
4.8600e-03 
7.0883e-03 
-1.8196e-01 
8.3131le-01 
-1.7869e+00 
-3.8659e+00 
4.3850e+00 
-5.5200e+00 
-7.6674e+00 
-1.9641le+01 
2.9900e+01 
-8.3855e+01 
-1.466le+02 
-8.7322e+00 


<3 -O0P73e-02 
-1.8258e-02 
-2.0696e-01 
9.6831le-01 
1.4301e+00 
-4.9158e-03 
4.3004e-01 
-7.2878e-02 
1.0061e-01 
-1.6629e-02 
3.6413e-02 
5 .8488e-03 
1.1409e-02 
-1.6150e-03 
-1.0547e-03 
-1.7207e-04 
7.9756e-03 
7.4508e-02 
-4.9034e-01 
-2.0083e-01 
2.-0019e+00 
-1.4247e+00 
1.6696e+00 
-2.9695e+00 
9.3531e+00 
-7.8291e+00 
-8.5556e+00 
1.1159e+00 
-1.3177e+02 


begf = 


-1.0822e+01 
6.8583e-01 
-4.7724e+01 
3.6012e+02 
4.7115e+02 
-6.4500e+01 
7.1903e+01 
-6.4391le+01 
6.4582e+00 
-1.7025e+01 
1.8515e+00 
1.2428e+00 
3.1364e+00 
-1.3148e+00 
-1.1379e-01 
-2.1182e-01 
8.7699e-01 
-7.4043e-01 
-2.6335e-01 
2.7014e-01 
-6.9940e-01 
2.4769e-02 
1.8978e-01 
-1.5582e-01 
2.3802e-01 
-4.7719e-02 
-3.1925e-02 
-3.5525e-03 
-3.9408e-03 


-2.7878e+00 
-4.1887e+00 
-2.8286e+01 
8.5027e+01 
1.4656e+02 
2.1242e+01 
6.8206e+01 
6.4952e+00 
1.9540e+01 
2.2886e+00 
7.2405e+00 
8.3604e-01 
1.3847e+00 
1.0517e-01 
-1.8892e-01 
-1.4546e-01 
6.0226e-01 
-5.0848e-01 
-1.8086e-01 
1.8552e-01 
-4.8031le-01 
1.7010e-02 
1.3033e-01 
-1.0700e-01 
1.6346e-01 
-3.2771e-02 
-2.1924e-02 
-2.4397e-03 
-2.7063e-03 


141 


Ccgt — 
Columns 1 through 6 


0 0 0 0 O O 
0 0 0 0 0 0 


Columns 7 through 12 


0 0 0 
0 0 0 0 0 0 


© 
oO 
oO 


Columns 13 through 18 


0 0 O -8.4999e-01 -1.5173e-01 -5.4154e-02 
0 O O -1.2095e-01 1.8226e+00 1.8269e+00 


Columns 19 through 24 


-3.7507e-O2 -3.8748e-01 6.6203e-03 4.5526e-02 -5.1326e-02 -1°4025e—a8 
-~8.6488e-01 -8.4670e-02 8.7789e=-01 <=-2.71176-02 2.5437e-0O1 -1.4015e-01 


Columns 25 through 29 


-6.3435€303 1.5178e-02 -3.9105€-02 -2.51136-02 "=e 77 [eee 
2.9992e-01 -5.7504e-02 -6.9062e-03 9.1004e-03 "=4°76205e-03 


dcgf 
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Closed loop state space representation of the H, 


limited-performance X-29 (unbalanced plant used in calculations). 


acgflim = 


Columns 1 through 6 


-1.0402e-02 -5.8688e-05 -1.0935e-01 1.0846e-03 9.5729e-02 -2.1473e-02 
4.7514e-05 -1.0127e-02 -2.2341e-02 2.5275e-01 7.5480e-02 -6.2276e-02 
-1.0953e-01 1.7080e-02 -2.1014e+01 -2.1657e+01 2.9586e+01 -8.7074e+00 
1.1830e-01 4.6756e-03 5.4776e+O1l -4.9923e+01 2.0305e+01 2.4992e+01 
1.0112e-01 2.5760e-02 3.047le+O1 -9.6706e+01 -7.0361e+01 5.0851e+01 
-7.0719e-03 7.5703e-02 -2.474le+00 -1.8064e+00 -2.0296e+01 -1.8520e+01 
-1.1890e-02 -4.0375e-03 -4.3794e+00 9.4786e+00 1.5457e+01 3.0890e+00 
3.5848e-02 -9.2905e-03 1.3787e+O1l -2.9232e+01 -5.0189e+01 1.1314e+0O1 
7.8006e-03 1.608le-02 2.7372e+00 -6.6157e+00 -1.1414e+01 -7.6888e+00 
6.9880e-03 -1.4086e-02 2.7883e+00 -5.6558e+00 -8.4502e+00 8.3607e+00 
3.6963e-03 -2.3286e-03 1.4203e+00 -3.0586e+00 -4.8787e+00 1.8586e+00 
9.8035e-04 -2.6757e-03 3.9874e-01 -7.6819e-01 -1.1181e+00 1.4626e+00 
-9.9671le-04 -1.2236e-03 -3.6218e-01 8.6485e-01 1.4612e+00 33777 7e-Cl 
4.9526e-04 -6.0858e-04 1.933le-0O1 -4.0393e-01 -6.2616e-01 3.8745e-01 
9.9360e-05 -1.1607e-04 3.8715e-02 -8.1164e-02 -1.2611le-01 7.4870e-02 
0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 

0 0 0 0 0 0 
eeroog7e+03 —-2.5519e+03 -3.2170e+04 3.2203e+04 3.5087e+04 -1.3450e+04 
=5.6119e+03 2.5455e+03 -2.7388e+04 -3.1437e+04 1.6463e+04 1.9671e+03 
-1.9979e+05 9.0622e+04 -9.7503e+05 <-1.1191e+06 5.8610e+05 7.0029e+04 
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Columns 7 through 12 


-4.0275e-03 
-2.4278e-02 
-1.7350e+00 
8.9719e+00 
8.5586e+00 
-1.2731e+01 
-7.5102e+00 
4.9757e+01 
5.8087e+00 
8.8626e+00 
4.3397e+00 
9.4872e-01 
-1.3423e+00 
§.5179e-01 
1.1123e-01 
0 


oo o0o 0 0 000 


© 


-4.1178e+03 
1.9746e+03 
7.0297e+04 


-2.6915e-02 
5.1044e-02 
-9.6190e+00 
-1.8751l1e+01 
2.7161e+01 
1.4809e+00 
-1.1428e+01 
-1.1776e+02 
-8.6183e+01 
-6.3838e+01 
-2.3127e+01 
-7.3025e+00 
6.0171e+00 
-3.4240e+00 
-6.8505e-01 
0 


oo oO 0 0 0 0 000 


-7.6282e+02 
-1.3742e+04 
-4.8921e+05 


3.4546e-03 
2.2025e-02 
1.5835e+00 
-8.5052e+00 
-8.5873e+00 
1.1718e+01 
1.0646e+01 
4.7311le+01 
-1.2232e+01 
-6.7146e+00 
-4.3150e+00 
3.8876e-01 
2.5122e+00 
-3.5627e-01 
-7.5761e-02 


oOo oO 0 00 0000 0 


3.6829e+03 
-1.8474e+03 
-6.5768e+04 


5.3706e-04 
-2.0440e-02 
-1.1432e-01 
8.1165e+00 
3.9686e+00 
-7.8909e+00 
-1.2192e+01 
5.8958e+01 
2.1820e+01 
-3.6223e+01 
-3.3307e+01 
-1.1983e+01 
2.5206e+00 
-4.1849e+00 
-8.2576e-01 
0 


oo oO 0 00 0000 


-2.4198e+03 
2.7284e+03 
9.7131e+04 
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-3.6531le-03 
2.5551le-03 
-1.3381e+00 
-9.6598e-01 
4.1757e+00 
-1.2714e+00 
-2.8839e-01 
-2.1270e+01 
4.0286e-01 
2.1952e+01 
-1.7274e+01 
-1.1844e+01 
7.0285e+00 
-4.9408e+00 
-9.9156e-01 
0 


oo oO 0 00 0000 


-6.5185e+02 
-1.3132e+03 
-4.6750e+04 


1.0496e-03 
2.5566e-03 
4.3285e-01 
-1.0186e+00 
-1.9185e+00 
1.6320e+00 
1.8363e+00 
-5.3409e-02 
-3.2960e+00 
8.7173e+00 
6.0655e+00 
-6.8722e+01 
1.3970e+01 
-5.8083e+01 
-9.0062e+00 


oo oO oO 0 00 000 0 


5.9980e+02 
-3.8383e+0O1 
-1.3665e+03 


Columns 13 through 18 


8.3492e-04 
-1.6840e-03 
2.8897e-01 
6.5715e-01 
-7.0494e-01 
-1.4062e-01 
-5.8267e-01 
7.5075e+00 
1.0924e+00 
-6.3426e+00 
9.2730e+00 
-8.7804e-02 
-1.3039e+02 
5.3211le+01 
1.3336e+01 
0 


oOo oO 0 QO 0 000 0 


1.1078e+01 
4.3911e+02 
1.5632e+04 


-5.3511e-04 
-4.2833e-04 
-2.0766e-01 
1.7396e-01 
7.8492e-01 
~4.9365e-01 
-4.5585e-01 
-1.7408e+00 
8.0532e-01 
-1.2524e+00 
-3.9489e+00 
5.6209e+01 
1.4771e+01 
-4.0347e+01 
-1.5833e+01 


O0o0O00000 000 O 


-1.9609e+02 
-9.0974e+01 
-3.2387e+03 


-8.1760e-05 
-1.6786e-04 
-3.3253e-02 
6.7004e-02 
1.4255e-01 
-~1.1507e-01 
-1.2596e-01 
-5.7699e-02 
2.2536e-01 
-5.5502e-01 
-4.7442e-01 
9.7762e+00 
1.1042e+00 
-~1.4513e+01 
-1.4561e+02 


0000000000 OC 


-4.2799e+01 
-9.6295e-01 
-3.4281e+01 


-5.2019e+00 
-~5.0395e-01 
-2.740le+01 
2.9927e+01 
2.6073e+01 
1.2523e-01 
-3.0793e+00 
8.7351e+00 
2.3598e+00 
1.3920e+00 
8.6600e-01 
1.7759e-01 
~2.8039e-01 
1.0853e-01 
2.1926e-02 
-4.1810e-01 
5.4740e+00 


oo oO 0 0000 0 


3.5657e+00 
3.7297e+02 
1.3278e+04 
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-1.1241e+00 
2.-3744e+00 
-7.7266e+00 
4.7994e+00 
2.-2229e+00 
-9.1880e+00 
-1.5243e-01 
2-9843e+00 
-1.4722e+00 
2.0218e+00 
4.6904e-01 
3.6550e-01 
9.0399e-02 
9.7591le-02 
1.8874e-02 

9.9600e-01 
-3.4240e-01 
0 


OOoO0O00 0 0 


oO 


1.3365e+00 
1.3980e+02 
4.9769e+03 


ooo 000 ogo 0 0 0 00000 


-2.2690e-02 
2-5850e+00 


oOo oO 0 00 000 


-1.4790e+03 
0 
0 


Columns 19 through 24 


oo oO 0 0 0 0 0 0 0 00000 
ooo 0o 0 0 0 0 0 0 020 000 
oOo 000 00 0 0 0 @lerere 
oo o 0 0 0 0 0 0 0 0000 
Oo oO QO 0 000 © © © Oe 
oo 0 0o 0 0 0 0 0 00000 


Oo 
Oo 
Oo 


-1.2130e-01 <-1.9480e-02 -9.4926e-04 4.4270e-05 -6.7120e-05 1.4510e-06 
-1.3860e+00 -1.0580e+00 3.8980e-03 -1.1640e-02 -6.3970e-03 -2.5090e-04 


0 0 1.0000e+00 0 0 0 
0 0 0 1.0000e+00 0 0 
0 0 0 0 1.0000e+00 0 
0 0 0 0 0 1.0000e+00 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 0 0 0 0 
0 0 -1.1430e+02 0 O -2.5290e+00 
-1.4910e+03 0 O -1.1490e+02 0 0 
O -5.3020e+04 0 O -1.8160e+03 0 
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Columns 25 through 29 


0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 
-2.1620e-05 -3.5400e-06 0 0 0 
-5.3620e-04 -2.9120e-04 0 0 0 
0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 

0 0 0 0 0 
1.0000e+00 0 0 0 0 
0 1.0000e+00 0 0 0 

0 0 1.0000e+04 0 0 

0 0 0 1.0000e+04 0 

0 0 0 0 1.0000e+04 

0 O -2.6970e+02 0 0 
-2.5360e+00 0 O -2.7010e+02 0 
O -1.7900e+01 0 O -6.0530e+02 
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begflim = 


5.2019e+00 
5.0395e-01 
2.7401e+01 
-2.9927e+01 
-2.6073e+01 
-1.2523e-01 
3.0793e+00 
-8.7351e+00 
~2.3598e+00 
-1.3920e+00 
-8.6600e-01 
-1.7759e-01 
2.8039e-01 
-1.0853e-01 
-2.1926e-02 


oOo O0oO0 00 00000 


-3.5657e+00 
-3.7297e+02 
~1.3278e+04 


1.1241e+00 
-2.3744e+00 
7.7266e+00 
-4.7994e+00 
-2.2229e+00 
9.1880e+00 
1.5243e-01 
-2.9843e+00 
1.4722e+00 
-2.0218e+00 
-4.6904e-01 
-3.6550e-01 
-9.0399e-02 
-9.7591e-02 
-1.8874e-02 


oOo oO0 oO 0 00 000 


-1.3365e+00 
-1.3980e+02 
~4.9769e+03 
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ccgflim = 
Columns 1 through 12 


0 0 0 0 0 0 
O 0 0 0 0 0 


Columns 13 through 24 


0 0 0 - 0 0 
0 0 O 0 1 0 


Columns 25 through 29 


0 0 0 0 0 

0 0 0 0 0 
dcgflim = 

0 0 

0 0 
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4. 


or 
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